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INTRODUCTION 


The operating procedure of a self-streamlining wind tunnel has already 
been discussed ^ and is summarised in the flow diagram shown on Figure 1. 

The iterative nature of the streamlining process/ requiring numerous measure- 
ments and calculations coupled with the need for a rapid and continual 
exchange of data between wind tunnel and computer, makes mandatory the use of 
a computer. 

The current operation of the Transonic Self-Streamlining Wind Tunnel 
(TSWT) involves on-line data acquisition with automatic wall adjustment. A 
tunnel run consists of streamlining the walls from known starting contours in 
iterative steps and acquiring model data. Each run performs what is described 
as a streamlining cycle. The associated control software is presented here. 

Development of the on-line control system for TSWT continues. It is 
anticipated that modifications to the existing software package will consist 
of minor changes to some of the subroutines. 

The introduction of further subroutines for 

(a) Calculation of model forces from the 'wall data* (wall pressures and 
positions) • 

(b) Solution of mixed imaginary flowfields with a shock present. 

(c) Prediction of wall shapes, to minimise boundary interference on a 
three dimensional model in a two dimensional test section. 

(d) Assessment of residual boundary interferences on a three dimensional 
model in a two dimensional test section. 

is planned as they become available. The modular architecture of the software 
will allow these additions to be made easily. 

The software package has been developed for simple application to other 
self- streamlining wind tunnels. The modular architecture allows individual 
control system s\±)routines to be utilised with existing users software. Also 
the software package is written in a general manner to minimise new application 
modification . 

During TSWT development mmerous programmes have been written to check 
sections of the control software. A nmber of these remain in use to assist 
with TSWT operation as follows: 
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i) Set both walls to known contours together or individually, 

ii) Allow operator modification of known wall contours, 

iii) Display current position of both walls, 

iv) Display and/or load contents of any specified data file record. 

In addition programs have been written to command a Tektronix 4662 plotter 
to display model pressure distributions, flexible wall Mach number distributions 
and wall shapes. 
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CONTROL SYSTEM HARDWARE OUTLINE 


An overview of the current control system is shown in Figure 2. Macro- 
scopically/ the system involves the interaction of the tunnel operator with 
the wind tunnel and computer to generate the required test data. 

The on-line control system has two distinct functions: 

1) to streamline the flexible walls which includes assessment of 
residual boundary interferences. 

2) to acquire and reduce test data from a model. 

These functions are achieved by two control loops between wind tunnel and com- 
puter linked to processing software for data manipulation. The control loops 
are for wall shape control and pressure data acquisition (Scanivalve control) . 
Data acquisition is automatic since the Scanivalve is stepped by computer com- 
mands and analogue pressure data is fed direct to the computer. 

The control system hardware has evolved about these two control loops using 
both analogue and digital data transfer. Microscopically / the control system 
becomes complex as shown in Figure 3. 

The system hardware consists of the computer and its peripherals communi- 
cating with a series of control and signal conditioning sub-systems housed in 
a control cabinet by the wind tunnel. This cabinet is then connected to the 
test section wall position sensors (linear potentiometers) and wall jack stepper 
motors^. 

The hardware will perform four functions: 

1) wall movement. 

2) wall and model pressure measurements. 

3) wall position sensing. 

4) system monitoring (not yet existing), 

The wall movement function involves the loading of forty motor latch boards with 
direction information (stop, forward-go or reverse-go) . Then the sending of 
a *go* pulse to the Pulse Sequence Generator (PSG) , starts actual wall movement. 
The PSG generates control pulses at a fixed frequency to the forty Motor Drive 
Boards (MDBs) . The MDB, using signals from the motor latch board and the power 
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supply, controls the sequence of power pulses transmitted to the S-phase stepper 
motor powering each wall jack. After a pre-determined and variable time interval 
the PSG is switched off and wall movement ceases. The wall movement control 
sequence is then repeated until satisfactory wall contours are achieved. 

The wall and model pressure measurements involve the driving of the 
scanivalve system by step pulses, to a required sampling port. Then the analogue 
signals from the pressure transducers are sampled by the computer after a suit- 
able settling time (i.e. 50 milli-seconds) . 

The wall position sensing function involves the sampling of forty analogue 
signals from the linear potentiometers attached to each wall jack mechanism. 

These signals are continuously available at the computer peripherals, but they 
have been found to be susceptible to electronic interference when the jack motors 
are switched on. Hence the wall position is only sampled with the walls station- 
ary. 

The system monitor is a necessary part of a practical digital control 
system. When this hardware becomes available it will provide information on the 
status of ‘hardware components, to allow quick error diagnostics during tunnel 
operations. 

The conditioned analogue signals from the test section (within the range 
1- 5 volt) are fed via the control cabinet to a DEC ADll-K Module for 12-bit analo- 
gue to digital coversion. This module is combined with a DEC AMll-K Multiplexer 
to make 6 4. channels available for input signals. 

The digital signals, at O and 5 volts , are transmitted to the wind tunnel 
in a code described in Appendix A. Device selection is by means of a 'telephone 
exchange* called the "address decoder**. In fact, commands to the test section 
are sent bo all devices but only one device is enabled to read the information, 
by address decoder selection. The versatility of a digital control system is well 
known and the reduction in interface wiring compared with an analogue system is 
significant. The address decoder has the capability of addressing 64 devices. 

The control system hardware layout in Figure 3 has been simplified for 
clarity. In practice there are mamerous synchronisation signal paths between 
devices, to ensure correct operation sequencing and to prevent 'race* problems. 

The tunnel operator monitors the control system from a computer VDU consol 
and inputs test parameters. The consol allows the rea^^l time display of test 
section and model data which is stored on the computer mass storage device. This 
data can subsequently be drawn in graphical foann on an XY Plotter (Tekronix 4662) . 



Hard copy of summarised run data is printed in real time on a line printer. 

The control of a model wake traverse mechanism has been incorporated in the 
system, utilising the wall movement technique already described. 

The control system hardware has excess capability with 11 spare digital 
address slots and 15 spare analogue input channels. This spare capacity may 
be used for the control of tunnel Mach number and model attitude in the future. 


; '.I 
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CONTROL SOFTWARE 


Computer software has been developed for the on-line TSWT control system 
using a versatile modular architecture. Hence the program has been reduced 
to a collection of manageable subprograms which can be combined to control the 
wind tunnel and output real time results or provide more detailed off-line re- 
analysis of previously acquired data. 

An overview of the control software package is shown below. 


File Type & Name 


File Storage Name 


Function 


Main Program 
(TSWT) 


OFLEX 


i) Control and sequence sub- 
routine calls. 

ii) Read test parameters from 
the operator. 


Subroutine 1 
(DATA) 


Acquire pressure data from the 
wind tunnel. 


Subroutine 2 
(REDUCE) 


Read tunnel data from disc storage 
and reduce raw pressure data from 
the wind tunnel. 


Subroutine 3 
(WAS) 

Subroutine 4 
(STAR) 


OJUDD 


Perform wall setting calculations, 

Calculate local boundary layer 
displacement thickness and Mach 
number along each wall. 


Subroutine 5 
(SOME) 


Assess wall induced interferences 
at the model. 


Subroutine 6 
(FORCE) 


OWING 

ONPL 


Calulcation of model forces for NACA 
0012-64 and NPL 9510 sections 
respective ly . 


Subroutine 7 
(SET) 


i) Store run data on disc 

ii) Output data to the teinninal 
and/or the plotter. 


Subroutine 8 
(WALL) 


OADJ 


Move the walls to new contours 


This breakdown of the software into modules has been extremely useful, parti- 
cularly for storage, editing and debugging purposes. 

The software written in FORTRAN IV language is run on a DEC 11/34 with 
a DEC RT-11 V4 operating system. The software is linked to a system library 
and a FORTRAN library to access functions and system subroutines and a Real 
Time System Library (RTSL) to access peripheral control subroutines. The com- 
plete compiled and linked program requires over 100 blocks (25.6K words) of 
memory spare. 

Current 16-bit computer processors are only capable of addressing 32K 
words (64K bytes) of real memory space. But of this, only 22K words is avail- 
able for a user's program, depending on the size of the operating system. 
Therefore to run the TSWT control software on a 16-bit machine a technique of 
overlaying has to be used, so that only part of the software is stored in the 
real memory at any instant during execution. 

Each subroutine is a self contained program communicating with the main 
program via common data blocks. So in theory only one subroutine is required 
in the real memory at any one time for execution. In practice, the subroutines 
have been grouped together to minimise the number of overlays, thereby reducing 
the time required for overlaying. The overlaying structure of the control 
software is shown below 


( 7762 words 


I 9470 words 


This program structure is implemented at 'link* time during the program 
generation cycle as below 



SEGMENT 1 

SEGMENT 2 

OVERLAY 

Subroutine 1 

Subroutine 5 

REGION 

Subroutine 2 

Subroutine 6 


Subroutine 3 

Subroutine 7 


7762 words 

Subroutine 8 



6558 words 

ROOT 

Main Program + System Library 

SEGMENT 

FORTRAN Library 


RTSL Library 
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R LINK 

*OFLEX = OFLEX/FORLIB,RTSL/C 
*OAD , ODR , OJUDD/0/ : 1/C 
*OERR r OWPL , OUT , OADJ/0 : 1 

These commands generated a runable program called OFLEX. The control software 
memory requirement was reduced from 27k words to 17.3k words. 

At *run* time the program OFLEX requires four data files to exist on the 
computer mass storage device. Data file ADC. DAT receives the raw analogue-to- 
digital counts of the 'wall and model data* for each streamlining iteration; 

PAD. DAT provides and receives sets of wall contours and the associated external 
imaginary wall velocity distributions; NPL.DAT or WING.DAT receives Cps from 
the NPL 9510 and NACA 0012-64 models respectively, for each streamlining itera- 
tion; TWST.DAT holds all fixed tunnel data, i.e. jack positions, potentianeter 
calibrations, scaling and coupling factors, matrix coefficients for camber inter- 
ference assessment and boundary layer information. RUN. DAT holds run data, i.e. 
ambient temperature and pressure, run number, iteration record number, and 
number of model tappings. The data files ADC. DAT, PAD. DAT and NPL . DAT/WING . DAT 
each hold 50 records: Records 1 to 3 in PAD. DAT hold data on the three aerodynami- 
cal ly straight contours (i.e. position and boundary layer thickness) . Records 1 
and 2 of NPL.DAT and WING.DAT files are used to store model tapping X and Y co- 
ordinates relative to the model leading edge. Records 4 to 50 are available to 
store data from each streamlining iteration. Hence iteration record mambers range 
from 4 to 50. When the iteration record number equals 50, ADC. DAT, PAD. DAT and 
NPL. DAT/WING. DAT must be copied to a data bank since the original data is then over 
written by subsequent iterations, starting with record 4. The upper limit on the 
iteration record number has been chosen to keep the data files in manageable propor 
tions (i.e. 25.6k words maximum size). The total storage requirement of data files 
to run OFLEX in 61.7k words. 

The data file RUN. DAT must be loaded with current run data before the control 
software is activated. This operation is performed by running a program OSTART, 
the software of which, called RUN, is described in section 3.11. 

A complete listing of the control software is described in the following 
sections. Where possible standard FORTRAN has been used but peripheral control 
commands are peculiar to the DEC system used. These subroutine calls can be 
grouped into Analogue to Digital sampling commands (ADC and RTS) , and programmable 
clock commands (SETR and LWAIT) . In addition there are calls to the system lib- 
rary routines (IPEEK and IPOKE) for digital input and output operations. 
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An example of the minimal printout from OFLEX associated with a typical 
streamlining cycle, involving two wall adjustments, is shown on Figure 5. The 
walls were initially set to contours stored in record 30 of file PAD. DAT and 
the run finished with wall contours as stored in record 32. The model had no 
pressure tappings and therefore no wing performance data was presented in the 
print-out. Note that the print-out from OFLEX is sent to a line printer (logi- 
cal Unit 7) while operator information is sent to a VDU Consol (Logical Unit 5) . 
The operator information consits of prompts to indicate the stages reached in 
the streamlining cycle and error warnings. 

The versatility of the software has allowed simple generation of programs 
for particular tasks such as tunnel data re-analysis. Using the existing sub- 
routines as building blocks, a new program has been made up of a series of these 
subroutines linked to a new main program. For example, data re-analysis is 
achieved by running the program ORLEX. The main program REAN is a modification 
of TSWT with different subroutine calls and an extended print-out demanded 
as described in section 3.10. The program structure is very similar to that of 
the control software and is generated using the following link command with a 
memory requirement of 17.1k words. 

R LINK 

*0RLEX = OREF,FORLIB,RTSL/C 
*0AD , ODR , O JUDD/0 : 1/C 
*0ERR ,0DST,0UT/0 : 1/C 

where QREF is the file storage name of program REAN 

An example of the extended printout is shown on Figure 6 for a typical re- 
analysis of raw TSWT data for one iteration of run 389. 

Should any new analysis technique become available, then a new subroutine 
could replace or supplement the existing subroutines, and be incorporated in 
OFLEX and/or ORLEX by minor adjustments to the main program and the LINK commands. 
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3.1. Program TSWT 


The main control program listed on Figure 4.1 reads tests parameters 
from the tunnel operator and sequences subroutine calls 


Lines 

0002 

- 0009 



Define all common data blocks used 
thus:- 

Block 

1 


NJ 


Effective number of jacks per wall 
(real plus dummy) . 




MT 

= 

Number of model pressure tappings. 




NR 


Number of tunnel reference samples 
during a scanivalve scan. 




CLl 


Printout control value. 




B1 

= 

Prandtl-Glauert scaling factor = 




PRl 


For development only. 




AKl 

= 

Aerodynamic coupling factor. 




AK3 


Wall movement scaling factor. 




AN 

= 

Model angle of attack (degrees) . 




R3 

= 

Chord Reynold s number. 




PP2 

— 

Dynamic pressure allowing for compres- 
sibility. 




ITRN 

= 

Run number. 

Block 

2 

Array 

RN 


Bottom wall potentiometer outputs, volts 



II 

RS 


Top wall potentiometer outputs, volts. 



II 

W 


Bottom wall imaginary wall velocities. 



II 

X 

= 

Top wall imaginary wall velocities. 



II 

P 


Top wall real wall static pressures. 



II 

Cl 

= 

Bottom wall real wall static pressures. 



II 

RD 

= 

Developnent data. 
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Block 3 

Array 

D 

= 

Measuring point co-ordinates. 


M 

WTY 

= 

Top wall movement from the straight 
wall contours (Inches) . 


II 

WBY 

= 

Bottom wall movement from the 
straight wall contours. 


II 

WL 

= 

Effective panel length per wall 
tapping . 


II 

B 

= 

Model Cps. 


II 

PD 

= 

Transducer Cps- 

Block 4 

Array 

E 

= 

Calculated imaginary external top 
wall velocities. 


II 

H 

= 

Calculated imaginary external bottcm 
wall velocities. 


II 

Y 

= 

Predicted top wall movements 
required (Inches) . 


II 

G 

= 

Predicted bottom wall movements 
required (Inches) . 


It 

WI 

= 

Real top wall velocities squared 
(V/Uco)^ 


II 

XI 


Real bottom wall velocities squared 

Block 5 

Array 

u 

= 

Difference between real and imaginary 
top wall velocities (u/y^) 


II 

V 


Difference between real and imaginary 
bottom wall velocities (u/„ ) 

* Uoo 


II 

DS 

= 

•straight' wall local 6* values. 


II 

DET 

=r 

Local A6* values. 


II 

CS 

= 

Matrix coefficients for camber inter- 
ference assessment. 
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Block 7 


Array 


PC 


Potentiometer calibrations { Volts/Inch) 


Block 8 


Block 9 

Line OOlO 
Line 0012 

Line 0013 

Lines 0014 - 0015 
Lines 0016 - 0022 


Line 0023 

Line 0024 

Line 0025 

Line 0026 - 0048 

Lines 0049 - 0059 
Line 0053 

Line 0057 


TAB 

AMP 

IFN 

lAM 


Array IX 


Ambient Temperature (Deg. c) 

Ambient Pressure (Cm Hg) 

Iteration Record No: 

Automatic In-file selection trigger 
(0-off; 1-on) 

Pressure channel zeroes. 


Initialise Arrays X & DET. 

Set streamlining trigger (PRl) to zero (i.e. walls 
not streamlined) . 

Set automatic in-file selection trigger (lAM) to 
zero (i.e. Non automatic selection) . 

Define File 5 as data file RUN. DAT. 

Read run data from RUN. DAT and load ITRN with the 
run number, IFN with iteration record number , TAB 
with ambient temperature, AMP with ambient pressure 
and MT with number of model taps. 

Set extended printout trigger (CLl) to zero (i.e. 
minimal print-out required) . 

NR = Number of tunnel reference samples in a pressure 
scan. 

NJ = Number of measuring points per wall (See Figure 
7) . 

Input test parameter: 

Angle of attack (AN) 

Subroutine call sequence. 

If no model tappings (i.e. MT=0) do not call Subroutine 
FORCE . 

If walls streamlined (i.e. PRl=l) do not call Sub- 
routine WALL 
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Lines 0061 - 0064 When the Iteration Record number is 50 and jump 

out of streamlining loop. 

Line 0066 Increment the Iteration Record number 

Line 0067 If walls are un-steamlined (i.e. PR1=0) continue 

the streamlining cycle. 

Lines 00 69 - 0078 Load run data array TUN and load file RUN. DAT with 

new values of TUN. 


3.2. Subroutine DATA (Incorporating Subroutines GAD and STEP) 

This subroutine primarily acquires all the tunnel pressures either from 
the tunnel via an automatically stepped scanivalve or from disc storage. Some 
data is then reduced. 

The option to acquire old data stored on disc allows re-analysis of pre- 
vious runs. All analogue to digital conversions associated with the pressure 
data are handled by software in Subroutine GAD, Automatic control of the Scani- 
valve is achieved by software in Subroutine STEP. The listing on Figure 4.2 can 
be broken down thus 


Lines 0002 - 0003 
Lines 0004 - 00011 

Line 0012 
Line 0013 

Line 0014 
Line 0016 

Line 0017 

Lines 0018 - 0021 


Dimension of data blocks DA, DB, DC and IDATA. 

Define only common data blocks required in Subroutine 
DATA. 

Initialise the first row of 2D- array IDATA. 

CHD = Model chord = 4 inches (NACA 0012-64) 

= 6 inches (NPL 9510) 

NJl = number of wall jacks. 

NR = number of tunnel reference samples in 
a pressure scan. 

IPl = Approximate low value of scanivalve position 

encoder output at pressure scan start condition. 

Transducer calibration valves 
(cm Hg per A-D count) . 
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Lines 0025 

Lines 0034 
Lines 0038 

Line 0041 

Lines 0044 


Lines 0053 

Line 0062 

Lines 0064 
Lines 0068 

Lines 0074 

Lines 0081 

Lines 0094 
Line 0102 
Line 0103 
Line 0104 
Line 0106 


0032 


0037 

0040 


0051 


0061 


0067 

0073 


0080 


0093 


0101 


If a re-analysis (i.e. CL1=500) / open data file 
4, select an existing record (IR) and fill IDATA 
with raw transducer A-D counts from a previous TSWT 
run. 

Load pressure channel zeroes into IDATA 

Identify all passes of the data acqination cycle other 
than the first and print a prompt. 

Initiate DO loop to acquire tunnel pressures using 
the scanivalve system. 

If the first data acquisition cycle (i.e. X(30)=0) 
read the four pressure transducer zero signals 
and halt the software until a *CR' is received from 
the terminal - transmitted when the tunnel flow is 
stabilised by the operator. 

Wait for a change in the encoder output generated 
by a scanivalve step using Subroutine STEP. If 
there is a step error halt program. 

Sample four transducer outputs using Subroutine 
GAD and load array IDATA. 

I 

Load pressure channel zeroes into array IX 

Open a data file 4 and fill a selected record (lOT) 
with 192 raw A-D count values from the iteration/ 
stored in IDATA. 

Convert raw data into pressures relative to tunnel 

I 

reference pressure (cm Hg) , and load into array PD. 

Load arrays P and Q with wall jack centreline pressures 
with dummy measuring points included (See Figure 7)* 

Load arrays P and Q with mid-jack centreline pressures* 

RD (NR +1) = test section total pressure reading. 

RD (NR + 2) = Ambient pressure (cm Hg) . 

RD (NR + 3) = Ambient temperature (deg C) . 

R1 = Stagnation pressure (cm Hg) 
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Lines 0107 - 0123 


Lines 0124 - 0140 


Line 0141 


Calculate tunnel reference Mach numbers at regular 
intervals through the pressure scan, and look for 
excessive fluctuations (i.e. Am > .Ol) • 

Calculate selected tunnel pressures for possible 
calibration checks. 

Store average tunnel reference static pressure in Bl. 


Subroutine GAD 
Line OOOl 

Line 0002 

Line 0003 
Lines 0004 - 0006 

Lines 0007 - 0011 

Lines 0012 - 0020 

Line 0021 

Subroutine STEP 
Line 0002 

Lines 0003 - 0005 
Line 0006 


Subroutine label with data transfer of value N and 
array IDATA. 

Dimension of IDATA (Average A-D counts) and IDT 
(Raw A-D counts) arrays. 

NS = Number of samples per channel 

Time delay of 50 milli-seconds using the programmable 
clock (Subroutine SETR) . 

Take NS samples of four analogue input channels at 
IkHz and store digital results in IDT. 

Convert all raw A-D counts from packed integer format 
to real nu m ber format (Subroutine CVSWG) then take 
the average of the NS samples per channel and store 
in IDATA. 

Stop the programracible clock. 


POKE digital output register with the command number 
16384 (See Appendix A) . Note no scanivalve address 
is required with the present hardware configuration. 

Time delay of 20 milli-seconds using the progreunmable 
clock . 

POKE digital output register with a zero value to 
clear. 
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3.3. 


Subroutine REDUCE 


This subroutine performs data reduction operations on tunnel pressure 
information and acquires further tunnel data from data files held on the mass 
storage device. The software listed on Figure 4.3. is described thus: 


Lines 0004 - 0010 

Lines 0011 - 0015 

Lines 0016 - 0025 

Line 0026 

Lines 0027 - 0041 
Lines 0042 - 0054 

Lines 0055 - 0080 
Lines 0083 - 0092 

Lines 0094 - 0096 

Lines 0097 - 0098 
Lines 0099 - 0109 

Line olio 


Define only common data blocks required in Sub- 
routine REDUCE. 

Set variables, as previously described for sub- 
routine DATA. 

Calculate the average freestream Mach number of the 
run. 

PP2 = ql compressibility correction to tunnel 

q- 

Calculate chord Reynold's number (R3) . 

Convert tunnel pressure data into Cp using the tunnel 
reference pressure associated with each group of six 
readings. 

Load model pressure coefficients into array B from 
transducer channels 1 and 3 data. Note this soft- 
ware is model dependent and configured for use 
of the NPL 9510 section. 

Specif iy IWF corresponding to 'Straight Wall" base 
data for Mach number panels M < .725, .725 <M> .825 
or M > .825 for re-analysis case only. Otherwise 
a diommy data base is used (i.e. IWF=1) . 

Load arrays DA, DB and PC with data from File 2 
(TUN. DAT) . 

Define File 3 as data file PAD. DAT. 

Define wall contour record number automatically 
(i.e. IAM=1) or by operator input. 

Load array DC with 'wall data' from File 3. 
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Lines 0111 

- 0118 

Load arrays OTY, WBY and DS with data from array 
DA. 

Lines 0119 

- 0124 

Load arrays CS and D with data from array DB 

Line 0125 


Read coupling factor (AKl) and scaling factor (AK3) 
from File 2, record 10. 

Lines 0126 

- 0129 

Calculate effective wall lengths between raid-jack wall 
points and store results in array WL. 

Lines 0130 

- 0139 

Load arrays RS, RN, X, W with data from array DC 
(See section 3.1 for array descriptions). 

Lines 0140 

- 0145 

Convert the wall Y co-ordinates into movement from 
straight contours in inches. 

Lines 0146 

- 0147 

Close data files 

Line 0148 


Signify the completion of the data acquisition cycle 
by setting X(30) = 1. 


3.4. Subroutine WAS 

12 3 

The wall adjustment strategy ' ’ (WAS) manipulates the calculated imbalance 
between real and imaginary wall velocities to generate a wall movement which 
will give zero wall loading or vorticity. 

The strategy requires interpolation of real wall velocities 
at regular intervals along each wall. To ensure accurate interpolation at the 
wall ends using curve fitting to the wall velocity distribution, two straight 
dummy wall extensions are added to each flexible wall, 15.25 cm (6 inches) and 
22.86 cm (9 inches) in length respectively. The dummy wall measuring points 
introduced are points 1 and 2 upstream where the wall velocity is assumed free- 
stream and points 22, 23 and 24 downstream, which have a wall velocity equal to 
that of measuring point 21 or jack 19. Measuring point 1 is also assumed to be 
the origin of the wall boundary layer. Measuring points 3 to 21 are real and 
correspond to the position of wall jacks. The software representation of each 
flexible wall, in terms of measuring point is shown in Figure 7. The software 
listed on Figure 4.4. is described below 

Lines 0002 - 0006 Define only comnon data required in subroutine. 
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Line 0007 
Lines 0009 - 
Lines 0011 - 

Lines 0016 - 


Lines 0040 - 
Lines 0042 - 
Lines 0045 - 

Lines 0054 - 

Lines ■ 0069 

Line 0099 
Lines 0104 - 


OOlO 

0012 

0038 


0041 

0044 

0053 


Dimension working arrays. 

Message to tunnel operator. 

Equalising top and bottom wall coupling and scaling 
factors. 


Compute external velocities for next wall shape where 


Q1 

TEMPI & TEPl 


U & V 


E & H 


XI & WI 


compressible dynamic pressure 
incompressible top and botton wall 
Cps. 

inbalance between real and imaginary 
wall velocities (Top and bottom walls) . 
external velocity perturbations 
computed for the next top and bottom 
wall shapes (v/uoo) 

2 

real wall velocities squared iv/u^) 


Load Z with inter jack streamwise spacing co-ordinates . 


Set up DO loop for both walls. 


Load A & XB with D U or V data in sets of fqur values 
for top or bottom wall. 


0067 


- 0098 


0137 


Cubic spline fit to each set of data to obtain wall 
velocity imbalance between jacks. 

Summation of velocity induced by the vorticity distri- 
buted along a wall at each jacking point. The result 
is stored in arrays S for the top wall and T for the 
botton wall. 

End of wall velocity analysis repeated for each wall. 

Numerical integration of top and bottom wall normal 
velocity components to generate jack movement demands. 
Raw movements are stored in TT and R. Arrays Y & G 
store the effective jack movements after scaling. 
Likewise E and H (the external velocity perturbations) 
are modified by scaling in lines 0133 & 0134. In 
lines 0135 & 0136 values in arrays Y and G are con- 
verted to compressible values using Bl = 


I 
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3.5. 


Subroutine STAR 


The function of STAR is to calculate Mach number and boundary layer dis- 
placement thickness values for each wall jacking point necessary only for re- 
analysis information. The boundary layer calculations use a numerical solution 
of the Von Karman Momentum Integral equation for a turbulent boundary layer. 

The software listing is shown on Figure 4.5 


Line 0002 
Lines 0003 - 0007 
Lines 0008 - OOll 
Lines 0012 - 0013 
Lines 0014 - 0015 
Lines 0016 - 0017 
Line 0018 
Lines 0027 - 0145 


Lines 0044 - 0045 

Lines 0046 - 0051 

Lines 0053 - 0062 

Lines 0063 - 0068 
Lines 0069 - 0077 

Lines 0081 - 0087 
Lines 0088 - 0092 


Dimension of array DELTA for 6* storage. 

Define only required common data blocks. 

Calculation of isentropic flow relationships. 

3 

Calculation of air density (kg/ft ) 

Calculation of air temperature (K) . 

Calculation of air viscosity (Ib/ft^^^) . 

CLl alue selects extended printout when equal to 500 

Overall Do loop performs calculations for each jack 
on both walls using sets of three jacks, labelled 
O, 1 and 2. 

Load XI and X2 with jack 1 and 2 co-ordinates relative 
to jack O. 

Ensure that the correct tunnel reference pressure is 
used with each wall pressure. 

Load SPl, SP2 and SP3 with top or bottom wall velocity 
perturbations . 

Calculate wall velocities UO, Ul and U2 (ft/sec) 

Calculate local wall mach number and store in P 
(top wall) or Q (bottom wall) . 

For the second measuring point assume a turbulent 

boundary layer growth according to 5* = .0213x 

^428 

Calculate the velocity gradient at jack 1 and store 
value in D2. 
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Lines 0093 - 0094 


Lines 0096 “ 0099 
Lines OlOO - 0105 


Lines 0106 - 0111 
Lines 0114 - 0122 

Line 0123 
Line 0127 
Lines 0131 - 0135 
Lines 0136 0141 

Line 0142 
Lines 0146 - 0163 


Guess the size of 5 * and store as PI. 

Calculate components of M.I. equation. 

Check if the measuring point is a special case 
(points 3 or 24 along each wall) and calculate SI 
with appropriate velocity gradient (SI = rate of 
change of 6* with streamwise distance) . 

ic if 

Calculate 5 from d6 and store value in P4. 

Home-in on correct value of PI (If P4 > PI increase 
PI by 25 E-6) . 

Load array DELTA with 6* values (inches) . 

Calculate A6* from ’straight wall’ to contoured wall. 

Output results to a terminal if CL1=500. 

Isolation of special cases requiring P4, ST & D2 
results to be stored in new addresses. 

Load array DELTA with 6* values (inches) . 

Calculate wall Mach numbers between jacks and store 
results in arrays P and Q/ corresponding to top and 
bottom wall data. 


3.6. . Subroutine SUME 

SOME assesses the residual interferences, due to any residual wall loading, 
along the tunnel centre-line and the model chord line. The three components of 
the interferences are related to estimated changes in Cj^. The quality of stream- 
lining is then assessed from these interferences and a measure of the average 
pressure imbalance across the walls. A listing is shown on Figure 4.6. 

Lines 0002 - 0007 Define only required Common Data blocks. 

Lines 0008 Define model in use IWT = O for NACA 0012-64 

IWT = 1 for NPL 951P 

Line OOlO CHD = Model chord (inches) 

Line 0013 OR = X co-ordinate of model pivot point which is at 

% chord and nominally mid way between the straight 
walls (Inches) . 
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Lines 0014 

Lines 0017 

Line 0020 

Line 0021 
Line 0023 

Line 0025 
Lines 0030 

Lines 0038 

Line 0046 
Line 0050 

Line 0051 
Lines 0054 

Line 0055 
Line 0057 


0016 


0018 


0037 


0043 


0074 


Calculate model chord line X and Y increments and 
store in AX and AH respectively. 

IPP = Position of model pivot point as a number of 
^ chords. 

OR = X co-ordinate of the model h chord point - the 
origin. 

Subroutine print trigger (OT) set to zero for print-out. 

Eliminate dummy jacks by reducing number of jacks 
per wall by four. 

Y3 = nominal tunnel semi-height (Inches) . 

Calculate imaginary wall velocity squared and sum 
the absolute value of the load difference between real 
and imaginary velocities squared for each jack. 

Store the results in EE (top wall) and F (bottom wall) . 

Calculate the average Cp difference between real and 
imaginary flows for each wall (called E) and store 
as ET and CEB, and print-out. 

If minimal print-out is specified (i.e. CLl = O) jump 
to line 0083 

For 37 points along the tunnel centreline summate 
the effect of wall voriticity, at the model and along 
the tunnel centreline. 

Xl = X co-ordinate of measuring point. 

Sum the velocity components of the wall vorticity 
at co-ordinate Xl on the tunnel centreline and store 
results in arrays UT and VT. 

X2 - Horizontal displacement between jack and analysis 
point (Inches) . 

Y2 = Displacement of bottom wall downwards from the 
centreline (Inches) . 
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Line 0065 


Lines 0067 


Lines 0076 
Lines 0083 

Line 0091 


Calculate the horizontal induced velocity perturba- 
tion component due to vorticity at measuring points 
n. (See Figure 7 ) 


U1 


Bottom 


n 

2tt 


Top 


2 2 


where Fn = Vortex strength 

= Local imbalance of wall real and imaginary 
velocities, 

= horizontal separation of analysis point 
and a vortex assumed at a measuring point 

n. 


Yn 

and 


vertical separation of same. 

panel length = distance between mid- jack 
points spanning measuring point n. 


In software notation 


U1 = — ^ *((U * R2) + (V * R3)) 
where WL, U and V are arrays. 

Calculate the vertical component of induced velocity 
perturbation VI due to vorticity at measuring points 
n 


0079 

0119 


VI 



Bottom 

* E — * r 

„ x^+y^ ^ 

Top n n 


If or = O print- out results. 

Calculate velocity perturbations for nine equidistant 
points along the model chord line. 

Calculate induced flow angles along the model chord 
line (AA2) . 
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Line 0092 
Line 0093 
Line 0103 - 0105 

Line 0107 
Lines 0110 - 0112 

Lines 0114 - 0118 

Line 0120 
Line 0123 
Line 0124 
Line 0128 
Line 0129 
Lines 0130 - 0139 

Lines 0140 - 0142 

Lines 0145 - 0154 

Line 0156 

Line 0160 - 0163 


Calculate induced Cp (CPI) , 

Store h chord point in address CP. 

Load Al and A2 with leading edge and trailing edge 
induced flow angles. 

SP = Contribution to induced camber. 

Differentiate between odd and even numbered measuring 
points . 

Load accumulators for Simpsons rule numerical integra- 
tion technique. 

CPE = Average Cp errors induced along model chord. 

P3 Induced camber in terms of Ac^^. 

A3 = Induced camber (degrees) . 

CL = Induced Angle of Attack in terms of AC^^. 

Al = Induced Angle of attack (degrees) . 

If a re-analysis print-out all residual interference 
information (i.e. CLl = 500). 

Print abbreviated version of residual interference 
information if CLl = O. 

Increment ISC for each quality of streamlining satis- 
fied. 

If all the qualities of streamlining are satisfied 
(i.e. ISC = 5) PRl = 1, the walls are streamlined. 

If the walls are streamlined print messages to the 
operator and the line printer. 


3.7. Subroutine FORCE (Based on a NASA numerical technique) 

FORCE integrates the model pressure distribution to give coefficients 
Cp and Cm* 
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3.7.1. 


NACA 0012-64 version - listed on Figure 4.7a. 


Lines 0002 
Lines 0007 
Line 0010 

Lines 0011 
Line 0013 

Line 0014 
Line 0015 
Lines 0016 

Lines 0024 

Lines 0031 

Line 0042 
Line 0043 
Line 0044 
Line 0045 
Lines 0047 

Lines 0051 

Line 0055 
Lines 0057 

Lines 0062 


0006 


0012 


0023 


0029 


0040 


0050 


0056 


0061 


0067 


Define required Common Data blocks. 

Dimension working arrays. 

WT = Number of model tappings +2 for dummy leading 
edge tappings. 

Open data file 2 (WING.DAT) 

Load array AG with X co-ordinates of model tappings 
(upper surface followed by lower surface) . 

Load AJ with Y co-ordinates as for the X co-ordinates. 

IHT = Number of model tappings per surface. 

Titling and test parameter input from the terminal 
(AN = angle of attack - degrees) . 

Calculate an effective vertical component for each 
model tapping, and store in array AH. 

Initiate a DO loop with special processing for first 
and last tappings on each surface and different anlysis 
for each surface. 

WF = Weighting Factor, in this case for Cj^. 

L6 = Local components of % for upper surface. 

L5 = Sum of Cn components. 

Store negative L6 in array A. 

Calculate the upper surface Cj^ component for the dummy 
leading edge tap. 

Calculate the component for the last downstream • 
tap on the upper surface. 

Cn “ Upper surface Cfj- 

Calculate the lower surface component for the dummy 
leading edge tap. 

Calculate % component for the last downstream tap on 
the lower surface. 
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Line 0066 

Lines 0068 - 0072 
Lines 0074 ~ 0088 
Line 0077 
Line 0078 
Line 0079 
Line 0083 
Line 0085 
Line 0089 

Lines 0098 - 0122 

Line 0103 

Line 0104 

Line 0105 

Line 0115 

Line 0116 

Line 0123 

Lines 0124 - 0131 

Lines 0132 - 0136 

Line 0137 - 0140 
Line 0141 


CNL = Lower surface 
L5 = Total model Cj^. 

Calculate Cj^ components for lower surface tappings. 

Do loop to calculate Cl components for each tapping. 

C5 = Upper surface C^ component. 

Store C5 in array AH. 

C3 = Sum of C^ components. 

CCP = Suuper surface C^. 

C5 = Lower surface C^ component. 

St = Lowers surface model C^. 

C3 = Total C . 

c 

Do loop to calculate C^ components about the model 
leading edge. 

BB = Upper surface Cj^ components. 

C9 = Sum of upper surface ^m components. 

TIT = Tapping position (% chord) . 

BB = Lower surface components. 

C8 = Sum of lower surface Cjjj components. 

STl = Total Cjjj. 

Output forces Cj^, C^, C^^^. 

Calculate coefficients Cl and C^ from Cj^/ C^ and 
angle of attack. 

Output wing performance coefficients. 

Close data file. 


3.7.2. NPL 9510 version - listed on Figure 4.7b 

This version of FORCE is the same as the NACA 0012-64 version in all respects / 
except for the number of model tappings (50 in total, with a split 32:18 on upper 
and lower model surfaces respectively) and the labelling associated with the print* 
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out. Also notice that print-out control commands ahve been inserted at lines 
0015 , 0092, 0108, 0112, 0120 and 0126. Setting CLl = O suppresses all print- 
out except data associated with the label *WING PERFORMANCE*. NPL 9510 model 
Cps are stored in NPL. DAT for subsequent plotting (See next section) . 


3.8. Subroutine SET 


SET controls the output of tunnel data to the VDU consol and the line printer 
and the output of results to disc storage. This software is configured for test- 
ing the NPL 9510 section. A listing is shown on Figure 4.8. 


Line 0002 
Line 0003 - 0009 
Line 0010 
Lines 0011 - 0012 

Lines 0013 - 0014 

Line 0015 
Line 0016 

Line 0017 
Lines 0019 - 0035 


Lines 0036 - 0039 
Lines 0042 - 0043 
Lines 0044 - 0053 


Lines 0055 - 0080 


Dimension data array DA. 

Define required Common Data blocks. 

Initialise some array elements for print-out purposes. 

Open wall output file 3 with 50 records each 512 
words in length (PAD. DAT) . 

Open wing output file 2 with 50 records each 128 words 
in length (NPL. DAT) . 

NJl = Number of wall jacks. 

Define model used: IWT = O for NACA 0012-64; 

IWT = 1 for NPL 9510 . 

CLl = O for minimal print-out, so jump to line 0055 

Output of pressure transducer Cps, external wall 
velocity for each jack and the predicted jack move- 
ment for zero wall interference. 

Print-out top wall labels. 

Define File 6 as WALM.DAT 

Output wall Mach number distributions stored in 
arrays P and 0 into an operator set (IMR) record of 
File 6 for subsequent plotting. 

Store potentiometer values for new wall contours in 
array DA. If CLl = 500 print-out information on jack 
X & Y co-ordinates, local Mach number, current poten- 
tianeter volt values, predicted potentiometer volt 
values for the next iteration. 
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Line 0060 


Lines 0061 - 0069 


Line 0078 

Lines 0081 - 0115 
Lines 0105 - 0114 


Line 0116 
Lines 0117 - 0118 

Lines 0122 - 0126 
Line 0128 
Lines 0129 - 0132 

Lines 0133 - 0134 


WM = Effective position of top wall contour allowing 
for boundary layer growth along the wall. 

If movement demands for the downstream jacks (No. 

16-20) exceed mechanical limits in terms of pot volts, 
set to safe minimum values. 

If CLl = 500 print out the local wall Mach Number bet- 
ween jacks. 

Repeat analysis and print-out for the bottom wall. 

Calculation of wall Mach number standard deviation 
for top and bottom walls for use in special re-analysis 
of tests with an empty test section. 

Set KT to equal the total mamber of jacks. 

Store top wall external velocity perturbations in 
array DA. 

Print prompts on wall and model output records. 

Write array DA into record IFN of File 3 on disc. 

Define the number of model tappings and store model 
Cps in record IFN of File 2 on disc. 

Close data files 2 and 3. 


3.9. Subroutine WALL (incorporating Subroutines INIT, START, MOVE and DIO) 

This subroutine controls the test section wall adjustments. Both walls 
are moved simultaneously in variable increments of movement. Each jack is 
commanded to move away from the model in its last adjustment to ensure the walls 
are rigid to air pressure loads. Numerous saftey checks are included to guard 
against jacks jamming or jacks out of control. 

Subroutine INIT initialises the control system. Subroutine START switches 
the jacks on and off for a time interval proportional to the average position 
error of the wall adjacent to the model pressure surface. Subroutine MOVE loads 
each jack with direction information and checks the data is loaded. Subroutine 
DIO handles all the Digital input and output operations between wind tunnel and 
computer. 

The software is dependant on the tunnel hardware and is listed on Figure 4.9 
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and described thus 


Line 0002 
Line 0003 
Lines 0004 


Line 0006 
Lines 0007 

Line 0011 

Line 0012 
Line 0013 
Lines 0014 

Line 0017 
Line 0018 - 

Lines 0019 * 

Line 0022 

Line 0023 
Lines 0024 • 


- 0005 


- 0010 


0016 


70 


0021 


0027 


Define Common data required. 

NWJ = total number of wall jacks. 

Define arrays thus 

PC = potentiometer calibration (volts per inch) 

IM = stores jack movement demands for checks. 

ADC = stores wall pot. volts for movement checks. 

PV = wall contour pot. volts. 

PL = minumum safe pot. volt readings for downstream 
jacks . 

Load array PL with values 

Define File 3 as PAD. DAT and load array PV with a set 
of wall contour pot. volts from record IFN of File 3. 

IF = 1 for the first pass through the wall adjust- 

ment cycle 

IP = anti-backlash overshoot increment. 

Initialise the control system. 

Load variables IW = Wall number (O : top/1 : bottom) 

ICNT = No. of attempts to move one jack. 

IS = No. of jacks correctly positioned. 

ITOL = Wall setting tolerance band in pot. volts. 

Calculate movement required for each jack from its 
current location and load each jack with direction 
information. 

Read the potentiometer output for jack (L) , convert 
the result to pot. volts, and store in IPV. 

IMOVE = Actual jack movement from previous position. 
Ignored on first pass. 

Store current jack pot. volts (IPV) in array ADC. 

If IPV > 980 (i.e. jack approaching mechanical limit) 
then stop. 
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Lines 0029 


Line 0038 


Line 0039 


Line 0041 


Line 0043 


Lines 0045 


Line 0051 


Line 0053 


Line 0055 


Line 0057 


0036 Single out downstream jacks (i.e. No. 16-20 and 36- 

40) and ensure that the demanded pot. volts are greater 
or equal to the minimum safety values . 

IPD = Demanded jack pot. volts adjusted to overshoot 
the wall position, towards the model, by an increment 
IP. 

IMV = Required jack movement in pot. volts. (+ve 
towards the model and -ve away from the model) . 

If the wall is being adjusted from its overshoot posi- 
tion (i.e. IP = O) only allow the wall to move away 
from the model (i.e. IMV < O) . Set IMV = O if the jack 
tries to move towards the model. This technique 
ensures the wall is rigid to air pressure loads by 
eliminating jack mechanism backlash in the correct 
sense. 

Change the sign of all bottom wall movement demands 
for hardware compatability. 

0049 Change the sign of top wall movement demands for all 

odd numbered jacks to accommodate the hardware con- 
. figuration for the top wall control. 

Jump to line 0063 if this is the first pass (i.e. 

IF = 1) . 

For jacks 1 and 21 do not check for jack jamming due 
to the small movement demands on these jacks positioned 
close to the wall anchor point (See figure 7) . 

If the demanded jack movement (IML(L)) is more than 
50 pot. volts and the jack potentiometer reading shows 
only a change of 10 or less pot. volts then warn the 
operator and abort the wall adjustment. 

If only one jack has yet to be correctly adjusted 
(IS = 39) increment ICNT for each pass through the 
adjustment cycle. 
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Line 0059 


Line 0061 


Line 0063 
Line 0064 

Line 0067 - 0068 
Line 0069 

Line 0071 
Line 0072 
Line 0074 
Line 0075 
Lines 0077 - 0080 


Line 0081 
Line 0082 - 0092 


Subroutine INIT 

Line 0002 
Line 0003 


If IS = 39 and ICNT = 6 then warn the operator and 
abort the wall adjustment. One jack may have moved 
out of its potentiometer measuring range. 

If the change in jack pot. volts between position 
samples is greater than 50 warn the operator and 
abort the wall adjustment. A jack may be moved the 
wrong way. 

Load array IM with values of required jack movement. 

If the required movement is less than or equal to 
ITOL go to statement 40 (line 0067) . 

Increment IS by one and set IMV = O. 

Load jack (L) with direction information using Subroutine 
MOVE. 

Set IF = O to indicate the first pass is complete. 

If IS = 40 the wall is correctly adjusted. 

Move the walls using Subroutine START. 

If IS less than 40 repeat the wall adjustment cycle 

If the wall has been correctly adjusted to the demanded 
contour (IP = O) finish. If the wall has been correctly 
adjsuted to the overshoot contour (IP = 5) repeat the 
wall adjustment cycle with IP = O. 

Set IW = O for top wall condition. 

For all downstream jacks (Nos. 16-20 and 36-40) warn 
the operator when their pot. limit has been reached. 


NWJ = 40 = Number of wall jacks. 

Set data operation as write before read IDI = 128. 
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Line 0004 - 0010 


Line 0007 


Subroutine START 
Line 0001 


Line 0004 - 0006 


Line 0007 


Line 0008 


Line OOlO 


Line 0012 - 0014 


Line 0016 - 0019 


Line 0020 - 0022 


Line 0023 


Subroutine MOVE 


Line 0001 


Line 0002 


Generate command numbers ICOM for a mother stop directive 
to each jack and send them individually to the wind tunnel 
using Subroutine DIO. 

Check the data has been correctly loaded. 


Subroutine statement with data transfer of the demanded 
jack pot. volts (PV) and the required movements or 
position errors (IM) . 

Determine the average position error of the bottom wall 
in pot. volts and store the result in lAM. 

RCT = number of clock counts between motor power on and 
off. 

The upper limit of RCT is 6CX). Set by wall damage 
considerations . 

The lower limit on RCT is 25. Set by resolution of the 
wall position measurements and speed considerations. 

If RCT > 50 load jacks 1 and 21 with a motor stop 
because they will overshoot by too large an amount. 

Generate a command number ICOM to start the motors 
(See Appendix A) and send it to the wind tunnel using 
Subroutine DiO. 

Variable time delay set by number of clock counts 
(RCT) at loo Hz. 

Stop all motors 


Subroutine statement with data transfer of IJ (Jack 
number) and IMV (required jack movement) . 

Determine whether IMV is negative, positive or zero 
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Lines 0003 - 0007 Load the jack direction data into IDR corresponding 

to either stop, forward-go or reverse-go. 

Line 0008 All data operations are write before read (IDI = 128) . 

Line 0009 Generate a command number ICOM. 

Line OOlO Send ICOM to jack IJ using Subroutine DIO. 


Subroutine DIO 

This subroutine performs the handshaking between computer and wind tunnel 
to transfer data on the digital I/O lines. 


Line 0001 

Line 0002 

Line 0003 
Line 0004 


Line 0006 
Line 0007 
Line 0009 

Line 0011 

Line 0012 


Subroutine statement with data transfer ICOM (Command 
number) ,IDI (Data operation) and INPUT (Data check) 

Load Digital output buffer (address 167774) with 
command number 45 to ensure all motors are switched off 

Send command number ICOM to the wind tunnel. 

Halt program until an external data accept signal is 
received on the digital I/O status register (hand- 
shaking between wind tunnel & computer) . 

Reset the Digital I/O status register to zero. 

If IDl = O (Write only operation) finish. 

Halt the program until an external data ready pulse 
is received on the digital I/O status register. 

Store the contents of the Digital input register in 
INPUT . 

Initialise the status register. 


3 . 10 Program REAN 


The main re-analysis program listed on Figure 4.10 reads test parameters 
and sequences subroutine calls. The program is very similar to TSWT, the main 
differences are: 
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1 - 


The printout trigger is set to give an extended print-out (i.e. 

CLl = 500) . 

2. All test parameters and run details are provided by the operator. 

3. The program contains no loops and will run through only once per 
activation. 

4. Subroutine WALL is not included. 

5. Subroutine STAR is included to provide additional information during 
the re- analysis. 


3.11 Program RUN 

This program has to be activated before each session of TSWT runs, to ensure 
the data file RUN. DAT contains the current ambient conditions, run and file 
numbers. 

The software is shown on Figure 4.11 and discussed below. 


Line 0002 
Lines CXD03 - 0004 
Line 0005 

Lines 0007 - 0009 

Lines 0010 - 0011 

Line 0012 
Line 0014 - 0021 

Lines 0022 - 0023 
Line 0025 
Line 0026 

Lines 0027 - 0028 


Define array TUN to hold run data. 

Define File 2 as data file RUN. DAT. 

Load array TUN with run data stored in RUN. DAT, record 

1 . 

Indicate if only an increment of the record number is 
required (Yes: INC = 1) . 

Read thermocouple A-D channel and load TUN (3) with 
the current ambient temperature. 

If INC = 1 jimp to line 0025. 

Input ambient pressure (cm Hg) , the run number and the 
file number as requested. 

Load TUN (1) & TUN (2) with new values . 

Increment the record number by one. 

Write the modified contents of array TUN to RUN. DAT 
on disc. 

Print the contents of array TUN as a check. 
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SYMBOLS 


wall velocity perturbation (V = U - U^o) - 
local wall velocity, 
freestream velocity. 

boundary layer displacement thickness. 

horizontal distance from the boundary layer origin (Inches) . 
Reynolds number based on boundary layer length (x) . 
vertical distance from the tunnel centreline (positive upwards) 
lift coefficient, 
pressure drag coefficient. 

pitching moment coefficient about the airfoil leading edge, 
normal force coefficient, 
chordwise force coefficient. 

Freestream Mach number . 

horizontal distance from leading edge (Inches) . 
vertical distance from the chord line (Inches) . 
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FIG. 1 SELF-STREAAALINING OPERATING PROCEDURE 
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FIG. 2 TSWT CONTROL SYSTEM OUTLINE 
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F5.gure 4 Listings of TSOT Control Software 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 

0019 

0020 
0021 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


SUBROUTINE WAS 

COMMON/ONE/NJfMTfNRfCLlfBl r PRl f AKl f AKS rANf R3»PP2 f ITRN 
C0MM0N/TW0/RN(30) fRS(30) fW(30) fX<30> fP(30) fQ(30) fRD< 20) 
C0MM0N/THREE/D<60) fWTY<30 ) fWBY<30) fWL(30) fB( 50) fPD(48r4) 
C0MM0N/F0UR/E(30) fH(30) fY(30) fG(30) fWI(30) fXK30) 
C0MM0N/FIUE/U(30) f V<30) f DS(60) fDET(60) fCS(20> 

DIMENSION A(4) fXB(4) fC<30f4) fZ( 30) fS(30> fT<30) 

REAL N 
WRITE(7f5) 

5 F0RMAT<///5Xf ' WAS COMPUTING NOW M') 

AK2 = AKl 
AK4 = AK3 
NJl = NJ-4 
36 MM = 1 

M = NJ 

DO 100 I = IfNJ 
IF(MM*EQ.NR> GO TO 10 
NT = (I-5)/6 
NT = NT * 6 

IF< (1-5) *EQ»NT) MM = MM +1 
10 01 = (RD(MM) - RD(NR + 1))/PP2 

TEMP = -P(I)/01 
TEMPI = Bl^TEMP 
TEMP = SORTCl-TEMPl) 

TEMP = TEMP-1 
U(I> = TEMP-X(I) 

Ed) = (AK3 jIcU(I)/2)+X(I) 

TEP = -Q(I)/Q1 
TEPl = B1*TEP ■ 

TEP = SQRT(l-TEPl) 

TEP = TEP-1 

0(1) = W(I)-TEP 

Hd) = Ud>-(AK4*0d)/2) 

WId) = <TEP+1))K(TEP+1) 

Xld) = (TEMP+1)*<TEMP+1) 

100 CONTINUE 

L = M-2 

DO 110 I = IfL 

110 Zd) = (Dd)+Dd + l))/2 
DO 175 NN = 1 f2 
NC = NN-1 
65 1=0 

35 DO 120 J = 1f4 

KI = I+J 
A(J) = Dd+J) 

15 IF (NC.EQ^O) GO TO 25 

XB(J) = Vd+J) 

GO TO 120 

25 XB(J) = Ud+J) 

120 CONTINUE 

00 = (XB(3)-XB(2))/(A(3)-A(2)) 

01 = XB(2)-00*A(2) 

02 = 1/(A(4)-A(1) ) 

03 = (XB(4)-00:<cA(4)-01)/((A<4)-A<2)):|c(A(3)-A(4))) 
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0050 

0051 

0052 

0053 

0055 

0056 

0057 

0059 

0060 
0061 
0063 

006 ^ 

0066 

0067 

0069 

0070 

0071 
0073 

0075 

0076 

0077 

0078 

0079 
TSWT 


CALL REDUCE 
CALL WAS 
CALL SUME 

IF(MT*EQ.0) GO TO 100 
CALL FORCE 
100 CALL SET 

IF(PR1*EQ»1*) GO TO 200 
CALL WALL 
200 IAM=1 

IF(IFN.EQ»50) WRITE(7f400) 

AOO F0RMAT(//60(')!c')/' REMEMBER TO COPY DATA FILES BEFORE NEXT RUN' 
C /60('J!c')) 

IF<IFN^EQ*50) GO TO 250 
IFN=IFN+1 

IF(PR1»EQ40*) go to 500 
250 TUN(1)=ITRN+1 
TUN(2)=IFN 

IF<IFN*EQ*50> TUN(1)=ITRN 
IF(IFN*EQ*50) TUN(2)=4* 

510 CALL ASSI6N<^f 'RUNtDAT' fOr 'OLD' f » ) 

DEFINE FILE A (6?256»UfIJR) 

WRITE<4'1) (TUN(M)rN==lf6) 

CALL CL05E(4) 

END 
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0001 


SUBROUTINE DATA 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 


0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 

0025 

0026 
0027 
0023 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 
0033 

0040 

0041 

0042 
,0044 

0046 

0047 

0043 

0049 

0050 

0051 

0052 


C 

C AUTO SCAN 
C 

DIMENSION DA ( 1 00 ) y DB ( 1 00 ) r DC ( 1 00 ) f SO ( 4 ) 

DIMENSION IDATA(49f4) 

COMMON/ONE/NJ»MTfNRfCLlfBlfPRlf AKlfAK3rANrR3fPP2f ITRN 
C0MM0N/TW0/RN(30) rRS(30) f W(30) f X<30) f P(40) f 0(40) r RD(20) 
COMMON/THREE/D (60) rWTY(30) f UDY( 30 ) > WL (30) ^ B(50 ) r PD(48f 4 ) 
C0MM0N/F0UR/E(30) fH(30),Y(30) f G(30) f WI (30) r XI (30) 
C0MM0N/FIUE/U(30 ) tV( 30) fDS(60) rDET(60) »CS(20) 
COMMON/SEUEN/PC ( 50 ) 

COMMON/EIGHT/TABr AMPf IFNf lAM 
C0MM0N/NINE/IX(4) 

DATA IDATA(lfl) »IDATA(lf2)f IDATA(lf3) f IDATA(l»4)/0T0y0f0/ 
C 

C DATA ACQUISITION 

C 

CHD = 4* 

NJl = NJ-4 
M = NJ 
NR = 8 
IPl = 2050 
S0(1) = *0331668 
S0(2) = *0142486 
S0(3) = *0134165 
S0(4) = .0141155 
IP0=0 

37 IF(CL1*EQ*0) GO TO 36 

lJRITE(5r31) 

31 FORMAT(/'$ DATA INPUT FILE = ') 

CALL ASSIGN ( 4 f f-1 f "OLD ' r » ) 

DEFINE FILE 4 (50f512fUf UR) 

WRITE(7f32) IFN 

32 F0RMAT(40(")Jc')/' ITERATION RECORD NO* = " f I4/40( " :|c' ). ) 

READ(4^ IFN) ( ( IHATA < N f J ) r N=1 f 49 ) f J==1f4) 

CALL CL0SE(4) 

GO TO 44 

36 IDATA(1f1)=IX(1) 

IDATA(1f2)=IX(2) 

IDATA(1f3)=IX(3) 

IDATA(1f4)=IX(4) 

IF(CL1*EQ*0*AND*X(30) .EQ.D WRITE (5f 800) 

800 FORMAT(" SCANIUALUE G0"/13('=")) 

DO 25 N=1f49 
IF(N*GT*1) GO TO 40 
IF(X(30) *EQ*1* ) GO TO 25 
CALL GAD(NfIDATA) 

WRITE(5f15) 

15 FORMAT(/" TURN ON WIND') 

READ(5f35) ISTOP 
35 F0RMAT(I2) 

GO TO 25 

50 F0RMAT(I3) 
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0070 
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0080 
0031 
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0033 
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0100 
0101 
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0105 


40 CALL STEP 

IENC=CVSWG( IADC<20f 1 ) » 1 ) 
IP2=IABS(IENC-IP1> 

IPl = lENC 

IF(IP2.GT»200) GO TO 30 
WRITE(7r43) 

43 FORHAT</' STEP ERROR') 

GO TO 180 

30 CALL GAEKNf IDATA) , 

25 CONTINUE 

IX<l)=IEiATA(lfl) 

IX(2)=inATA(lf2) 

IX(3)=IDATA(1?3) 

IX<4)=IDATA(lf4) 

CALL ASSIGN<4r 'ADC*DAT' »0f 'OLD ' f » ) 
DEFINE FILE 4(50r512»U» IJR) 

WRITE<7r42) IFN 

42 FORMAT(/' ITERATION RECORD NO. = ' rIA) 

WRITE(4'IFN)((IDATA<Nf J)fN=lf49)r J=l74> 
CALL CLOSE (4) 

C 

C DATA REDUCTION 

C 

44 DO 45 J = 1»4 

DO 5 N = 2»49 

IDATA<N»J)= IDATA(1»J)-IDATA(N»J> 

NN = N-1 

PD(NN»J)= SO<J)»IDATA(Nr J) 

5 CONTINUE 

45 CONTINUE 

DO 110 N = 1»2 
P(N) = 0.0 
0(N) = 0.0 
110 CONTINUE 

DO 115 N = 22r24 
P(N) = PD<19»2) 

Q(N) = PD(19j4) 

115 CONTINUE 

DO 120 N = 3r21 

lA = N-2 

P(N) = PD<IA»2) 

Q(N) = PD(IA»4) 

120 CONTINUE 

IT=21 

DO 210 N=25t33 
IT=IT+1 

IF(N.EQ,28) IT=IT+1 
P(N)=PD(IT»2) 

Q(N)=PD(IT»4) 

210 CONTINUE 

RD<NR-H) = PD<48»1> 

RD(NR+2) = AMP 
RD(NR+3) = TAD 
ADI = 0. 
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0123 

0124 
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0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 
DATA 


R1 = RD(NR+2) - RD(NRil) 

DO 55 J = 1»NR 
ITR = 1 + ((J-1):JC6) 

RD(J) = PD(ITR»1) 

ADI = ADI + RD(J) 

R2 = RD(NR+2)-RD(J> 

R3 = 0*28571«AL0G(R1/R2) 

PPl = 5.0 t (EXP<R3)-1) 

TFM = SORT (PPl) 

IF(J.EQ.l) TFM1=TFM 
TR = RD(J)/2.54 
IF(J.EQ.l) GO TO 55 

IF(ABS(TFM-TFM1) .GT.O.Ol) WRITE<7»56) 
56 FORMAT(/" MACH NO. ERROR V20( ) 

55 CONTINUE 


HI 


PD(1j-2) 

+ 1 

RD(1) 

H2 


PD(43»2) 

+ 

RD(6) 

H3 

= 

PD (20 1.2) 

+ 

RD(4) 

H4 

S5 

PD(lr4) 

+ 1 
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H5 

ss 

PD(43i.4) 

+ 

RD(6) 

H6 

ss 

PD(20»4) 

•f 

RD(4) 

H7 

.•= 

PD(2»1) 



HI 


Hl/2.54 



H2 

=: 

H2/2.54 



H3 

sr 

H3/2.54 



H4 

= 

H4/2.54 



H5 


H5/2.54 



H6 


H6/2.54 



H7 

=2 

H7/2.54 



AS 

22 

ADI /NR 



TR 

22 

AS/2.54 



SP 

22 

RD(NR+l)/2 

♦ 54 


B1=AS 

180 RETURN 
END 
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000.1 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
GAD 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 
STEP 


SUBROUTINE GAD <N?IDATA)- - 

DIMENSION IDATA(49s-4)fIDT(500) 

NS=15 
TEND = 0 

CALL SETR(4r0f50«rIEND) 

CALL LWAIT<0f lEND) 

I END = 0 
IRTS==0 

CALL RTS<IDTt60j »NS^ 16>4» »2» IRTSr IDUM) 
CALL SETR(4?9)-l»rIEND) 

CALL LWAIT(Of IRTS) 

DO 5 J " lr4 
SUM = 0 

DO 70 NN = IrNS 

IS = ((NN-l)>!c4)+J 

IDT (IS) = CVSWG(IDT(IS)rl) 

SUM = SUM + IDT (IS) 

70 CONTINUE 

IDATA(MjJ) = SUM/NS 
5 CONTINUE 

CALL SETR(-l»y?) 

RETURN 

END 


SUBROUTINE STEP 

CALL IPOKE( •167774^16384) 

IEND=0 

CALL SETR(4»0r20. 7 IEND) 
CALL LWAIT(OrlEND) 

CALL IPOKE( " 167774 »0) 

RETURN 

END 
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0001 


SUBROUTINE REDUCE 

0002 


DIMENSION DA(100)VDB<100) rDC(lOO) »S0(4) 

0003 


DIMENSION IDATA(49»4) 

0004 


COMMON/ONE/NJfMT»NRfCLlrBl>PRlf AKlf AK3»ANfR3fPP2f ITRN 

0005 


C0MM0N/TU0/RN(30)»RS(30)fW<30)fX(30):>P<40)fQ(40)»RD(20) 

0006 


C0MM0N/THREE/D(60) rUTY<30) /WBY<30)VWL(30) f B(50) fPD(48r4) 

0007 


C0MM0N/F0UR/E(30) fH(30) f Y(30) f G(30) f WK30) f XI (30) 

0008 


C0MM0N/FIVE/U(30) >0(30) f DS(60) fDET(60) fCS<20) 

0009 


COMMON/SEVEN/PC (50) 

0010 


COMMON/EIGHT/TABfAMPfIFNfIAM 


u 

C 

DATA ACQUISITION 

0011 

u 

AS=B1 

0012 


CHD = 4* 

0013 


NJl = NJ-4 

0014 


M = NJ 

0015 


NR = 8 

0016 


MM = 1 

0017 


R1 = RD(NR+2)-RD(NR+l) 

0018 


R2 = RD(NR+2) - AS 

0019 


R3 == 0*2857lJ«cAL0G(Rl/R2) 

0020 


PPl = 5*0)«c(EXP(R3)-l) 

0021 


AMI = SQRT(PPl) 

0022 


B1 =SQRT(1-PP1) 

0023 


WRITE(5»90) AMI 

0024 


URITE(7r90) AMI 

0025 

90 

FORMAT(//' MACH NO^ = "fF8*4) 

0026 


PP2 = 1*0 + (0*25>KPP1) + (0»025)KPPl>lcPPl) 

0027 


CCl = 1 + (0^2>KPP1) 

0028 


CC2 = 2*5 « ALOG(CCl) 

0029 


CC3 EXP(CC2) 

0030 


D1 = RD(NR+2)){(8 4 998E-3/(2734l5+RD(NR+3)) 

0031 


D1 = D1/CC3 

0032 


AT = RD(NR+3)+273*15 

0033 


AT = AT/CCl 

0034 


CC4 = 0*92j{^AL0G( AT/273 *15) 

0035 


01 = 11*52:4^EXP(CC4) 

0036 


PP3 = (AS~RD(NR+1) )/PP2 

0037 


UO = SQRT(56*369>!cPP3/Dl) 

0038 


R2 = UO >!c D1 >!« 3218E4/V1 

0039 


R3 = R2 jK CHD/12* 

0040 


WRITE(7f95) R3 

0041 

95 

FORMAT(//" REYNOLDS NO* - 'fE12*6) 

0042 


DO 60 J = 2f4 

0043 


MM=1 

0044 


DO 65 N = lf48 

0045 


IF(MM*EQ*NR) GO TO 170 

0047 


NT = (N-5)/6 

0048 


NT = NT :4c 6 

0049 


IF((N“5)*EQ*NT) MM = MM + 1 

0051 

170 

Q1 = (RD(MM)“RD(NR+1) )/PP2 

0052 


PD(NfJ) = -PD(NfJ)/Q1 

0053 

65 

CONTINUE 
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0039 

0090 
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60 CONTINUE 

DO 500 J=16f24 

IJ=J+8 

JB=J+29 

IJB=J+22 

B(IJ)=PD(J,3) 

IF(J^GT*21) GO TO 500 
B(JB)=PD(IJBf3) 

500 CONTINUE 
MM=1 
NP=1 

DO 600 K=2f42 

NT=(K-5)/6 

NT=NTX«6 

NJJ=<K-7)/6 

NJJ=NJJ>i£6 

IF ( ( K-5 ) . EQ * NT ) MM=MM+1 
IF((K“7)tEQ*NJJ) GO TO 600 
Ql= ( RD ( MM ) ~RD < NR+ 1 ) ) /PP2 
B(NP)=(RD(MM)~PD<Kf 1) )/Ql 
NP=NP+1 

IF(K«EQ»28) NP=33 
600 CONTINUE 

CALL ASSIGN(2f 'TUN*DAT' »0f 'OLD' » r ) 

DEFINE FILE 2 ( 12 r256f U f UR) 

IF<CL1*EQ*0*) IUF=1 
IF(CL1»EQ*0*) GO TO 77 
URITE<7f75) 

75 FORMATC/' FILE NO* ?'/' 1 FOR M=<*725'/' 2 FOR *725<M>* 825' ) 

WRITE(7f76) 

76 FORMAT('$ 3 FOR M=>*825 ANS = ' ) 

READ(7?35) IWF 

35 F0RMAT(I2) 

77 REAri(2'IUF) <DA< J) r J=1 r80) 

READ(2'11) (DB<J) U=lr33) 

NJ2 = NJl + NJl 
READ(2'12) (PC(J)>J=lrNJ2) 

CALL ASSIGN(3r 'PAD*DAT' yOy 'OLD' > f ) 

DEFINE FILE 3 <50T512>Uf INR) 

IF(IAM*EQ*1) 60 TO 20 
WRITE(5yl40) 

140 FORMAT(//'$ WALL RECORD? = ') 

READ(5f35> IR 
IF(CL1*EQ*0*) GO TO 41 
GO TO 40 
20 IR=IFN-1 

41 URITE(7f30> IR 

30 FORMAT(/" WALL CONTOURS RECORD ='fI4> 

40 READ<3'IR) (DC( J) f J=1 r 83) 

DO 150 K = 1»80 
IF(K.LE*20) WTY<K> = DA(K) 

IF(K*LE*40.AND*K*GT*20) WBY<K-20) = DA(K) 

IF(K*GT*40) DS<K“40) = DA(K) 

150 CONTINUE 
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0119 


DO 155 J = If 33 

0120 


IF(J.LE.9) CS(J) = DB(J) 

0122 


IF(J.GT*9) D(J-9) = DB(J) 

0124 

155 

CONTINUE 

0125 


READ(2'10) AKlfAK3 

0126 


DO 156 J = IfNJl 

0127 


NJJ = J + 2 

0128 


WL(J> = <D(NJJ+1>-D(NJJ-1) )/2 

0129 

156 

CONTINUE 

0130 


DO 160 K = If 88 

0131 


IF(K»LE.20) RS(K) = DC(K) 

0133 


IF<K*GT*20*AND.K*LE*40) RN(K-20) 

0135 


IF ( K . GT « 40 ♦ AND » K . LE ♦ 64 ) X ( K-40 > 

0137 


IF(K.GT*64> W(K-64) = DC(K) 

0139 

160 

CONTINUE 

0140 


DO 175 J = IfNJl 

0141 


WTY( J)=UTY<J)-RS< J) 

0142 


WTY( J)=WTY( J)/PC( J) 

0143 . 


WBY< J)=RN<J)-WBY(J) 

0144 


WBY ( J ) =WBY ( J ) /PC < J+N J1 ) 

0145 

175 

CONTINUE 

0146 


CALL CL0SE(3) 

0147 


CALL CL0SE<2) 

0148 


X(30)=l, 

0149 


RETURN 

0150 


END 

REDUCE 



= DC(K) 
= DC(K) 
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0001 


SUBROUTINE WAS 

0002 


COMMON/ONE/N J r MT t NR f CLl f B1 r PR 1 f AK 1 r AK3 f AN r R3 f PP2 f I TRN 

0003 


C0MM0N/TW0/RN(30)fRS<30)fW(30) fX(30) fP(30) fQ<30) fRD<20) 

0004 


C0Mf10N/THREE/D<60) fWTY(30) fWBY(30) fWL(30) fB(50) fPD(48f4) 

0005 


C0MM0N/F0UR/E<30) fH(30) fY(30) fG<30) fWK30) fXK30) 

0006 


C0MM0N/FI0E/U(30) fM(30) fDS(60> fDET(60) fCS(20> 

0007 


DIMENSION A(4)fXB(4)fC(30f4)fZ(30)fS<30)fT(30> 

0008 


REAL N 

0009 


WRITE(7f5) 

0010 

5 

FORMAT ( ///5X F ' WAS COMPUTING NOW NO 

0011 


AK2 = AKl 

0012 


AK4 = AK3 

0013 


NJl = NJ-4 

0014 

36 

MM = 1 

0015 


M = NJ 

0016 


DO 100 I = IfNJ 

0017 


IF<MM*EQtNR) GO TO 10 

0019 


NT = (I“5)/6 

0020 


NT = NT * 6 

0021 


IF(<I-5)*EQ,NT) MM = MM +1 

0023 

10 

Q1 =; (RD(MM) - RD(NR + 1))/PP2 

0024 


TEMP = "P<I)/01 

0025 


TEMPI ~ Bl^TEMP 

0026 


TEMP = SQRT(l-TEMPl) 

0027 


TEMP = TEMP-1 

0028 


U(I) = TEMP-X<I) 

0029 


Ed) = (AK3>i<U(I)/2)+X(I) 

0030 


TEP = -Q<I)/Q1 

0031 


TEPl = Bl)!cTEP ' 

0032 


TEP == SQRT(l-TEPl) 

0033 


TEP == TEP-1 

0034 


U(I) = W(I)-TEP 

0035 


H<I) « W<I)-(AK4JlcO(I)/2) 

0036 


WId) = <TEP+1)5 |c(TEP+1) 

0037 


Xld) = (TEMP+l)>lc(TEMP+l) 

0038 

100 

CONTINUE 

0039 


L = M-2 

0040 


DO 110 I = IfL 

0041 

110 

Zd) = (Dd)+Dd + l))/2 

0042 


DO 175 NN = 1f2 

0043 


NC = NN-1 

0044 

65 

1 = 0 

0045 

35 

DO 120 J = 1f4 

0046 


KI = I+J 

0047 


A<J) = D<I+J) 

0048 

15 

IF (NC*EQ*0) GO TO 25 

0050 


XB<J) = Od+J) 

0051 


GO TO 120 

0052 

25 

XB<J) = Ud+J) 

0053 

120 

CONTINUE 

0054 


UO = <XB<3)-XB<2))/<A(3)-A(2)> 

0055 


VI = XB<2)-V0>lcA(2) 

0056 


V2 = 1/(A<4)-A(1)) 

0057 


V3 = (XB<4)-V0JKA(4)-Vl)/((A(4)-A<2))5*c(A<3)-A<4>>) 
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0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0067 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 
0090 
0031 
0082 

0083 

0084 

0085 

0086 

0037 

0038 

0039 

0090 

0091 

0092 140 

0093 

0095 

0096 

0097 45 

0098 130 

0099 175 

0100 75 

0101 


V4 = (XB(1)-V0*A(1)-V1)/(<A<1)-A(2))*<A(3)-A(1))) 
06 = 02* (03-04) 

05 = 04-06*A(l) 

1 = 1+1 

PI = A(2) + A(3) 

C(t»l) = 01-A(2)*A(3)*05 
C(If2) = 00+05*Pl-06*A(2)*A(3) 

C(If3) = 06*Pl-05 
C(If4) = -06 

IF (I*LT.(M-3)) GO TO 35 
LI = M-2 
DO 130 J = 2rLI 
ZO = Z(J) 

Z02 = Z0*Z0 
Z03 = Z02*Z0 
SS = 0 
K = M- 3 
DO 140 I = IfK 
Y1 =D(I+1) 

CO = C(Ifl) 

Cl = C(If2) 

C2 = C<Ip3) 

C3 =? C(lfA) 

Y2 = D(I+2) 

Y2SQ = Y2 * Y2 
YISQ = Y1 * Y1 

50 = C0+C1*Z0+C2*(Z02)+C3*(Z03) 

TEMP = ABS(Y2-Z0)/ADS<Y1-Z0) 

51 = ALOG(TEMP) 

52 = <C1+C2*Z0+C3*Z02)*(Y2-Y1) 

53 = (C2+C3*Z0)*((Y2S0)-(YlSQ))/2 

54 = C3*< <Y2SQ*Y2)-(Y1SQ*Y1) )/3 

55 = SS+S0*Sl+S2+33+S4 
CONTINUE 

IF (NCtEQ^l) GO TO 45 

S(J) = SS/6»28319 

GO TO 130 

T(J) = SS/6t28319 

CONTINUE 

CONTINUE 

R = 0 

TT = 0 


0102 


S(l) = 0.0 


0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
0111 
0112 
0113 


T(l> = 0.0 
DO 150 I =1>NJ1 
TO = S(I) 

RO = T(I) 

T1 = Z(I> 

11 = I+l 
T2 = D(I1> 

T3 = S(Il) 

R3 = T(I1) 

T4 = Z(I1> 

12 = 1+2 


Figure 4.4.2 



0114 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 
WAS 


T5 = D(I2) 

T6 = S<12) 

R6 = T(I2) 

T7 = Z(I2) 

FSl = (T6-T3)/(T7r*T4) 

FS2 = (R6-R3)/(T7-T4) 

T2SQ = T2)!cT2 
T5SQ = T5>!cT5 

T8 = (FS1-(T0-T3)/(T1-T4) )/(T7-Tl) 

T9 = FSl - T8 * T7 

P2 = ( T3-T9>KT4 ) « ( T5-T2 ) +T8* < ( T5SQ*T5 ) - ( T2SQ>lcT2 ) ) /3 
TT = TT+P2+(T9-T8)KT4)#((T5SQ)-(T2SQ))/2 

R8 = (FS2-(R0-R3)/(T1-T4))/(T7-T1) 

R9 = FS2 - R8 «T7 

P3 = ( R3-R9HCT4 ) t ( T5-T2 ) +R8>lc < < T5SQ)5cT5 ) - ( T2SQ)icT2 ) ) /3 

F = E(I2) 

R = R+P3+(R9-R8>!cT4))!c((T5SQ)~(T2SQ))/2 
Y(I2) = <AK35!(TT)KAK2)!{AK45|{R) 

G(I2) = (AK4>i<R) + (AKl!icAK3>KTT) 

E(I2) = E(I2) + ((H(I2)-W(I2)+V(I2))«AK2> 

H<12) = H(I2) + ((F-Ua2)-X(I2))!{cAKl) 

Y(T2) ~ B1 ^ Y(I2) 

8(12) = B1 >!c 8(12) 

150 CONTINUE 
RETURN 
END 


Figure 4.4.3 
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SUBROUTINE STAR 
DIMENSION DELTA(60) 

C0MM0N/0NE/NJfMTfNRFCLlfBlfPRl»AKl»AK3fANrR3,PP2fITRN 
C0MM0N/TU0/RN(30) rRS(30)fW(30) f X(30) rP<30) » Q(30) f RD<30) 
C0MM0N/THREE/D(60) rWTY(30) rWBY(30) rWL(30) fB<50) 

COMMON/FOUR/E (30) fH(30)fY(30) f G(30> f WI (30 ) fXI ( 30) 
C0MM0N/FIVE/U(30)fV(30)fDS(60) f DET(60) »CS(20) 

PPl = 1 - (BDIcBl) 

CCl = 1 + (*2*PP1) 

CC2 = 2.5 * ALOG(CCl) 

CC3 = EXP ( CC2 ) 

D1 = RD(NR+2)*8.998E“3/(273.15+RD(NR+3)> 

D1 = D1/CC3 

AT = RD(NR+3)+273.15 

AT = AT/CCl 

CC4=.92)!cAL0G(AT/273^15) 

U1 = 11.52>JfEXP(CC4) 

IF(CLl.EQ.O.) GO TO 56 
URITE(7f5) 

FORMAT (///5Xf" DELTA STAR CALCS.'//' TOP WALL'/) 

WRITE (7f 10) 

F0RMAT(2X>' TAP N0.'»2Xf' DU/DX'f5X»' MACH N0.'y6Xf< D*',9X' DD*') 

KJ I*? =5 N I * ^ 

RRl = RD(NR+2) - RD(NR+1) 

MM = 1 

DO 2 N = 2rNJ2 
NN = N 

IF(N.GT.NJ) NN = N-NJ 
IF(NN.GT.(NJ-2)) GO TO 2 

IF(CL1.EQ.500. .AND.N.EQ. (NJ+D) WRITE (7f 11) 

FORMAT(//' BOTTOM WALL'/) 

IF(CL1.EQ.500. .AND.N.EQ. (NJ+D) WRITE (7>10) 

IF(N.EQ.(NJ+D) MM = 1 
IF(N.EQ.(NJ+D) GO TO 2 . 

IF (N.E0.NJ.0R.N.EQ.NJ2) GO TO 12 
XI = D(NN) - D(NN-l) 

X2 = D(NN+1) - D(NN-l) 

IF(MM.EQ.NR) GO TO 120 
NTS = N-5 
NT = NTS/6 

NT = NT»6 

IF(NT.EQ.NTS) MM=MM+1 
IF(N.GT.NJ) GO TO 110 
SPl = P(N) 

SP2 = P(N+1) 

SP3 = P(N“1) 

GO TO 100 
NN = N-NJ 
SPl = Q(NN) 

SP2 = Q(NN+1) 

SP3 = Q(NN-l) 

PP3 = (SP1+RD(MM)-RD(NR+1))/PP2 
PP4 = (SP2+RD(MM)-RD(NR+1))/PP2 
PP5 = (SP3+RD(MM)-RD(NR+1) )/PP2 


Figure 4.5.1. 



0066 

0067 

0068 


U1 = SQRT<56»369*PP3/D1) 
U2 = SQRT<56«369)JcPP4/Dl) • 
UO = SQRT<56*369*PP5/D1) 


C 

C LOCAL MACH NO. CALCS 


0069 

C 

R2 = RD(NR+2) - (SP1+RD(MM>) 

0070 


R3 = 0. 28571 ?!cAL0G(RRl/R2> 

0071 


PPl = 5.0!jc(EXP(R3)-l) 

0072 


AMI = SORT (PPl) 

0073 


IF(N.GT.NJ) Q(NN-1)=AM1 

0075 


IF(N.GT.NJ) GO TO 23 

0077 


P(NN-1)=AM1 

0078 

23 

IF (N.EQ.2.0R.N.EQ.(NJ+2>) 60 TO 22 

0080 


GO TO 33 

0081 

22 

R2 = U05KDDK32.18E6/V1 

0082 


R1 = <R2>!CX1)/12 

0083 


Cl = 0.142857«AL0G(R1) 

0084 


C2 = EXP(Cl) 

0085 


C3 = 0.00127 * XI 

0086 


P6 = C3/C2 

0087 


GO TO 133 

0088 

33 

Y1 = U1 - UO 

0089 


Y2 = U2 ~ UO 

0090 


A1 = (Y2-((X2>!:Yl)/Xl))/(<X2*X2)-(Xl>!cX2)) 

0091 


B2 = (Yl-(Al){c(Xl)|cXl)))/Xl 

0092 


D2 = (2>!cAl>!(Xl) + B2 

0093 

12 

DO 3 Ml = If 1000 

0094 


PI = (25 + M1>!(100)/1E6 

0095 


NODE =1 

0096 

122 

P2 = (Ul)icPl)!<Dl3}:32.18E6)/Vl 

0097 


C4 = 0.25 « AL0G(P2) 

0098 


C5 = EXP(C4) 

0099 


P3 = 0.0128/C5 

0100 


IF (N.ECJ.NJ.0R.N.EQ.NJ2) GO TO 66 

0102 


SI = P3-(3.45*{PlM(D25icl2/Ul) 

0103 


IF (N.EQ.3.0R.N.EQ. (NJ+3)) GO TO 77 

0105 


GO TO 88 

0106 

66 

SI = P3-(3.4#Pl>!cD35!a2/Ul) 

0107 


P4 = P5+(0.5J!c(Sl+S2))icXl/12> 

0108 


GO TO 99 

0109 

77 

P4 = (S1#X1/12> + P6 

0110 


GO TO 99 

0111 

88 

P4 = P5+(0.5«(S1+S2)#X1/12) 

0112 

99 

IF (N0DE.EQ.2) GO TO 155 

0114 


IF (P4.LT.P1) GO TO 111 

0116 

3 

CONTINUE 

0117 

155 

IF (P4.GT.P1) GO TO 144 

0119 

111 

PI = PI - 25E-6 

0120 


NODE =2 

0121 


GO TO 122 

0122 

144 

PI = PI + 25E-6 

0123 


DELTA(N) = 16.8 *P1 . 

0124 


N2 = N-2 


Figure 4.5.2. 



0125 

0127 

0128 
0129 

0131 

0132 

0133 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 


0146 

0147 

0148 

0149 

0151 

0152 

0153 

0154 

0156 

0157 
0159 
0161 
0162 
0163 

0165 

0166 
STAR 


IF(N.GT*NJ) N2 = N-6 
DET(N2) = DELTA(N) ~ DS(N2) 

NNJ = NN-2 

IF(CL1»EQ»0») GO TO 36 

URITE(5»25) NNJrD2»AMl»DELTA(N) rDET(N2) 

25 F0RMAT(5Xf I2r4F12»4) 

IF(N«EQ*<NJ2-2)) WRITE(5f35) REY»FPD 

35 FORMAT (/' UNIT REYNOLDS NO* = D>lc FPG = ' fF8*A/> 

36 P5 = P4 
32 = SI 

IF (J-(NJ-1)> 2»44»55 
55 IF (J~(NJ2-D) 2»44»2 

44 D3 = D2 

GO TO 2 

133 DELTA(N) P6 * 16*8 
REY = R2/12 
FPD = DELTA(N) 

2 CONTINUE 

C 

C WALL MACH NOS BETWEEN JACKS 
C 

IWALL=0 

350 DO 300 J:=25r33 

SF*1 ^F* ( J ) 

IF<IWALL*NE*0) SP1=Q(J) 

R2=RD(NR+2)-(SPl+RD(5)> 

R3=0 ♦ 28571>!(AL06 ( RR1/R2) 

PPl=5*0Jlc(EXP(R3)-l) 

IF(PP1*LE*0) GO TO 300 

AM1=S0RT(PP1> • 

IF(IWALL.EQ*0) P(J)=AM1 
IF(IWALL*NE*0) Q(J)=AM1 
300 CONTINUE 

IWALL=IWALL+1 
IF<IWALL.EQ*1) GO TO 350 
RETURN 
END 


Figure 4.5.3. 
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SUBROUTINE SUME 

COMMON/ONE/N J f MT ^ NR f CLl r B1 f PRl f AK 1 f AK3 f AN f R3 » PP2 f ITRN 
C0MM0N/TU0/RN<30) fRS(30) f W(30) fX(30) fP(40) fQ( 40) fRD<20) 
COMMON/THREE/D ( 60 ) f UTY ( 30 > f WBY ( 30 > f WL ( 30 ) f B < 50 ) f PD ( 48 f 4 ) 
COMMON/FOUR/E (30) fH(30 ) f Y(30> fG<30) f WI (30 > fXI (30) 
C0MM0N/FIUE/U(30) f0(30) fDS(60) fDET(60> fCS(20) 

DIMENSION UT(37)fVT(37) 

IWT=1 

TOPI = 6*283185 
CHD = 4. 

IF(IWT*EQ*1) CHD=6*0 

0R=24*56 

ANl = AN)|c. 01745 

AX = (CHD*C0S(ANl))/8*O 

AH = (CHD^SIN(ANl))/8.0 

IPP = 7 

IF(IWT*EQ*1) IPP = 4*33 
OR = OR - (dPP-3)#AX) 

OT = 0 
TOT =0*0 
NJl = NJ-4 
Y4 = 0*0 
Y3 =3*0 
MN = 0*0 
EE = 0* 

F = 0* 

IJ = NJ-2 

DO 160 I = 3fIJ 

01 = W(I)+1 

V2 = X(I)+1 

W(I) = VI >lc VI 

X(I) = V2 « V2 

EE = EE + ABS(XI(I)-X(I)> 

F = F + ABS(WI(I)~W(I)) 

160 CONTINUE 

GET = EE/NJl 
CEB = F/NJl 

IF<CL1*EQ*0* )URITE(5f170) CETfCEB * 

WRITE(7f170) CETfCEB 

170 F0RMAT(/10Xf ' WALL CP ERROR V5Xf' TOP - 'fF8*4f 

C 5Xf' bottom - 'fF8*4) 

IF(CL1*EQ*0*) GO TO 95 
WRITE(7f5) 

5 F0RMAT(10Xf' residual ERRORS') 

WRITE(7f10) 

10 F0RMAT(7Xf' X'fSXf' U/UFS'f6Xf' V/UFS') 

DO 161 I = 1f37 

XI = OR + ((I-19)*1*0) 

56 SUMU = 0*0 
SUMV = 0*0 
DO 162 J = IfNJI 
X2 = XI - (D(J+2)) 

Y1 = Y4+Y3-WTY(J)-DET(J) 

Y2 = Y4-Y3-WBY(J)+DET(JTNJ1) 


Figure 4.6.1 



0058 


PI = X2 3|c X2 

0059 


P2 = PI + (Y1*Y1) 

0060 


P3 = PI + (Y2*Y2) 

0061 


R1 = X2/P2 

0062 


R2 = Y1/P2 

0063 


R3 = Y2/P3 

0064 


U1 = (U<J+2)*R2)+(0(J+2)»R3) 

0065 


U1 = (U1)KWL<J))/T0PI 

0066 


01 = (U(J+2)+0(J+2))*Rl 

0067 


01 = <01)|cWL(J))/T0PI 

0068 


SUMU = SUMU + U1 

0069 


SUMO = SUMO + 01 

0070 

162 

CONTINUE 

0071 


IF (MN^GT*0) GO TO 98 

0073 


UT(I) = SUMU 

0074 


OT(I) = SUMO 

0075 


X3 = (1-19) t 1*0 

0076 


IF (0T*GT*0> GO TO 161 

0078 


WRITE(7fl5) X3 fUT(I)>0T(I) 

0079 

15 

F0RMAT(3F12*4) 

0080 

161 

CONTINUE 

0081 


WRITE <7^20) 

0082 

20 

F0RMAT(/5Xf' model ERRORS' f25X 

0083 

95 

SO = 0*0 

0084 


SE = 0^0 

0085 


SM = 0*0 

0086 


DO 164 JJ = lf9 

0087 


XI = OR - (<3“JJ)*AX) 

0088 


Y4 = -(IPP-JJ):«tAH 

0089 


MN = 1 

0090 


GO TO 56 

0091 

98 

AA2 = ATAN(SUM0/(1+SUMU>) 

0092 


CPI = 1*0 ~ (d+SUMU):«t<l+SUMU) 

0093 


IF (JJ*EQ*3) CP = CPI 

0095 


TOT = TOT + CPI 

0096 


IF <0T*GT*0) 60 TO 104 

0098 


XI = XI " OR 

0099 


IF<CLl*EQ»Od GO TO 104 

0101 


WRITE(7f25) XlrSUMU^SUMOfCPl 

0102 

25 

F0RMAT<4F12.4) 

0103 

104 

IF (JJ.EQd) A1 = AA2 

0105 


IF (JJ*EQ*9) A2 = AA2 

0107 


SP =: AA2 Jjc (l-CS(JJ)) 

0108 


IF (JJ*EQ*l40R*JJ*EQ*9) 60 TO : 

0110 


MB = JJ/2 

0111 


MB = MBJJC2 

0112 


IF (MB^EQ^JJ) GO TO 101 

0114 


SO = SO + SP 

0115 


GO TO 164 

0116 

101 

SE = SE + SP 

0117 


GO TO 164 

0118 

102 

SM = SM + SP 

0119 

164 

CONTINUE 

0120 


CPE = TOT/9^0 


Figure 4.6.2 



0121 

0122 

0123 

0124 

0125 

0127 

0128 

0129 

0130 

0132 

0133 

0134 

0135 
0136. 

0137 

0138 

0139 

0140 
0142 

0144 

0145 

0146 
0148 
0150 
0152 
0154 
0156 
0158 
0160 
0162 

0163 

0164 

0165 
SUME 


PI = SM+(2*S0)+(4*SE) 

P2 = (A)</3))«cPl. 

P3 = 2tP2 

A3 = (A1-A2) * 59^29578 
IF(CLl.NE^O) WRITE(7>35) 

35 F0RMAT(/10Xf ' EFFECT' f25Xf ' DELTA CL') 

CL = TOPI * A1 
A1 = A1 >!c 59*29578 
IF(CL1*EQ*0) GO TO 58 
WRITE(7f30) AlfCL 

30 F0RMAT(/5Xf' ALPHA ERROR = 'fF8*4f' DEGREES' fSXfFS *4) 
WRITE(7f40) A3fP3 

40 F0RMAT(/5Xf' INDUCED CAMBER = '»F8*4f' DEGREES' f2XfF8*4) 
WRITE<7r45) CP 

45 F0RMATC/5Xf' UEL*ERROR CP =-‘'»F8*4) 

WRITE (7^55) CPEf-CP 

55 F0RMAT(5Xf ' AVERAGE' » 6Xr '=' rF8*4f 13X» F8*4///> 

58 IF<CL1*EQ*0) URITE(5f 57) AlfA3rCPE 

IF(CL1*EQ*0) URITE(7f57) AlfA3fCPE 
57 FORMAT(' RESIDUALS = 'f3F8*4) 

ISC=0 

IF(CET*LE* *01)ISC=ISC+1 
IF<CEB*LE* *01)ISC=ISC+1 
IF(ABS(A1)*LE* *015)ISC=ISC+1 
IF<ABS(A3) *LE**07)ISC=ISC+1 
IF(ABS(CPE) ♦LE**007)ISC=ISC+1 
IF<ISC*EQ*5)PR1=1* 

IF<CLl*EQ*0*AND*PRl*Ea*l* ) URITE(5r800) 

IF(PR1*EQ*1*) WRITE<7f810) 

800 FORMAT(/6('5K')f 'TURN OFF WIND' f 6< '*' ) ) 

810 F0RMAT(//9X>23< '#' )/9X7 'ft WALLS STREAMLINED ft'/9Xf23< 'ft' > ) 

RETURN 
END ’ 
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SUBROUTINE FORCE 

C0MM0N/0NE/NJfMTfNRfCLlfBlfPRlFAKlfAK3fANfR3fPP2FlTRN 
C0MM0N/TU0/RN(30)fRS(30 ) fW( 30) fX(30)FP<30) r Q<30) »RD(20) 
C0MM0N/THREE/D(60) »WTY(30) rWBY<30 ) fWL(30) fB<50) f PD( 48f4) 
COMMON/FOUR/E (30 ) fH(30) fY(30> fG(30) fWI (30) fXI(30) 
C0MM0N/FIVE/U(30) fV(30) fDS(AO) fDET(60) fCS(20) 

DIMENSION A(50)fBB(50)fAJ(50)fAG(50)fAH(50) 

REAL L5fL6fL7fJ 
INTEGER S 
MT MT+2 

CALL ASSIGN(2f "WING.DAT' fOf 'OLD' f f > 

DEFINE FILE 2 ( IOOf 128fUf INR) 

READ(2'1) (AG(L)fL=1fMT) 

READ(2'2) (AJ(L)fL=1fMT) 

IHT = MT/2 
WRITE(7f30> 

F0RMAT(///15Xf' NACA SECTION ANALYSIS'/20Xf ' 0012-64'/ 

C /'$ RUN NO. = ') 

URITE(7f7) ITRN 
WRITE(7f9) . 

FORMAT(/'^ ALPHA = ') 

WRITE(7f12) an 
F0RMAT('+'fF6.2) 

F0RMAT('+'fI5) 

DO 54 K = 2f(MT-1) 

AH(K) = (AJ(K+l)-AJ(K-l))/2 
AH(1) == AJ(2)/2 
AH(IHT+1) = AJ(IHT+2)/2 
AH(IHT) = (.0012 - AJ(IHT-l))/2 
AH(MT) « (.0012 - AJ(MT-l))/2 
L5 =0 

DO 3 IQ = IfMT 
IF (IQ.EQ.l) GO TO 540 
IF (IQ.EQ.IHT) GO TO 580 
IF (IQ.EQ.(IHT+1)) GO TO 630 
IF (IQ.EQ.MT) GO TO 530 
IF (IQ.GT.(IHT+1)) GO TO 780 
WF = (AG(IQ+l)-AG(IQ-l))/2 
L6 = UFtBilQ) 

L5= L5-L6 
A(IQ) = -L6 
GO TO 3 

WF = AG(IQ+l)/2 
L5 = -UF«B(IQ) 

A(IQ) = L5 . 

GO TO 3 

WF = (1.0 - AG(IHT-l))/2 
L6=WF>«cB(IQ) 

L5=L5-L6 
A(IQ)=-L6 
CN = L5 
GO TO 3 

WF = AG(IQ+l)/2 
L6 = WFJKBdQ) 


Figure 4.7a,l 
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0080 
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0101 
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L5=L5+L6 
A<IQ) = L6 
GO TO 3 

530 WF = (1*0 - AG<MT-l))/2 
L6 = WF*B(MT) 

L5 =L5+L6 
A(MT)= L6 
CNL = L5-CN 
GO TO 3 

780 UF=( AG( IQ+1 )-A6( IQ-1 ) >/2 

L6 = UF#B(IQ) 

L5 = L5+L6 
A(IQ)= L6 
3 CONTINUE 
890 C3 = 0 

no 40 lY = IfMT 
IF <IY*GT*(IHT+1)> GO TO 980 
C5= (B(IY)*AH(IY)> 

• AH(IY)=C5 
C3 = C3+C5 

IF <IY*EQ*IHT) GO TO 970 
GO TO 40 
970 CCP = C3 
GO TO 40 

980 C5«(B<IY)*AH(IY)) 

AH<IY> = C5 
C3 = C3 + C5 
40 CONTINUE 

ST = C3-CCP 

IF<CLl*EQ.O*) GO TO 196 
WRITE (7fl65) 

165 FORMAT(////flOXf ' UPPER SURFACE') 

WRITE <7fl94) 

194 F0RMAT(2Xf' %CH0RD'f3Xf' CP LOCAL' f3Xf' CN LOCAL' f3Xf 
C ' CC LOCAL' f3Xf' CM LOCAL') 

196 C9=0^ 

C8=0* 

DO 5 S=lfMT 

IF (S*EQ^(IHT+D) GO TO 1240 
IF (S^GT*<IHT+1)> GO TO 1270 
BB(S) = (-A(S)#AG(S)) + (AH(S)3|cAJ(S)) 

C9 = C9+BB(S) 

TIT = AG(S)5!cl00 

IF(CLl*EQ.O*) GO TO 5 

WRITE (7f 185) TIT>B(S)fA(S)fAH(S)fBB(S) 

GO TO 5 

1240 IF<CL1*EQ*0*) GO TO 1270 

WRITE<7rl95) 

195 F0RMAT(///10Xf ' LOWER SURFACE') 

WRITE (7fl94) 

1270 BB(S) = (A(S)5KAG<S) > + ( AH(S):(cAJ(S) ) 

C8= C8-BBCS) 

TIT = AG<S)3!c 100 
IFCCLltEQ^O*) GO TO 5 
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FORCE 


WRITE <7fl85) TITFB(S)fA(S)fAH(S)f-BB(S) 

1 85 FORMAT < 4X » F4 ♦ 1 » 4F1 2 ♦ 4 ) 

5 CONTINUE 

STl = C8+C9 
URITE(7f205) 

205 F0RMAT(//18Xf 'PRESSURE' »2Xf 'SUCTION' f2Xf 'TOTAL'/) 

WRITE (7^206)CNfCNLfL5 

206 F0RMAT<12Xf'CN'f2Xf3F9*4) 

WRITE(7f207)CCPfSTfC3 

207 F0RMAT(12Xr 'CC'f2Xf3F9*4) 

WRITE<7f208)C9fC8rSTl 

208 F0RMAT<12X>'CM'f2Xf3F9#4) 

AT = AN 0.017453 

CSl = COS(AT) 

SN = SIN(AT) 

CL = (L5*CSl)-(C3)fcSN) 

CD = (C35ICCS1) + <L5)kSN) 

WRITE (7f225) 

225 F0RMAT(///22Xf ' WING PERF0RMANCE'/18Xf 'CL' r lOXf 'CD' f IOX f 'CM' ) 
WRITE (7f235) CLfCDfSTI 

235 FORMAT <10Xf3F12. 4) 

CALL CL0SE<2) 

END 
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SUBROUTINE FORCE 

COMMON/ONE/NJfMTfNRfCLlfBlfPRlf AKl»AK3f AN>R3»PP2f ITRN 
C0MM0N/TW0/RN(30)fRS<30)fU(30)fX(30)fP(30)fa(30>fRD<20) 
C0MM0N/THREE/EK60) fWTY(30) f UBY <30 ) > WL(30 ) fB( 50) fPD(48f4) 
COMMON/FOUR/E < 30) »H(30) r Y(30) rG(30) fWI( 30) fXK30) 
C0MM0N/FIVE/U(30) fV(30) r EiS(60) f DET(60) tCS<20) 

DIMENSION A(50)fBB(50)rAJ(50)rAG(50)>AH<50) 

REAL L5»L6fL7fJ 
INTEGER S 

CALL ASSIGN(2f 'NPL^DAT'rOr 'OLD'f f) 

DEFINE FILE 2 <50f 128fUf INR) 

READ(2'1) (AG(L) fL=1fMT) 

READ<2'2) (AJ(L)fL=1fMT) 

IHT = 32 

IFCCLl^EQ^O*) GO TO 31 
WRITE<7 f30) 

F0RMAT(///15Xf' NPL section analysis V20Xf ' 9510 V 

C /'$ RUN NO* = '> 

URITEC7 f7> ITRN 
WRITE<7f9) 

FORMATC/'^ ALPHA = ') 

WRITE<7f12) an 
F0RMAT('+'fF6*2) 

F0RMAT('+'fI5) 

B<31)=(B(30)+B<32) )/2* 

DO 54 K = 2f(MT-1) 

AH<K) = <AJ(K-M)-AJ<K-1))/2 
AH(1) = AJ<2)/2 
AH(IHT+1) = AJ(IHT+2>/2 
AH(IHT) = (*0024) 

AH(MT) = (*0024) 

L5 =0 

DO 3 IQ = IrMT 
IF (IQ*EQ*1) GO TO 540 
IF (IQ*EQ*IHT) GO TO 580 
IF (IQ*EQ*(IHT+1>) GO TO 630 
IF (IQ*EQ*MT) GO TO 530 
IF (IQ*GT*<IHTfl)) GO TO 780 
IJF = (AG(IQ+l)-AG(IQ-l))/2 
L6 = WF>KB(IQ) 

L5= L5-L6 
A(IQ) = -L6 
GO TO 3 

WF = AG<IQ+l)/2 
L5 = -UFJKBdQ) 

A(IQ) = L5 
GO TO 3 

WF = ♦5-<AG(IHT-l)/2) 

L6=UF:«(B(IQ> 

L5=L5-L6 
A(IQ)=-L6 
CN = L5 
GO TO 3 

WF = AG(IQ+l)/2 
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L6 = UF*B(IQ) 

L5=L5+L6 
A(IQ) = L6 
GO TO 3 

530 WF = ♦5-(AG(MT-l)/2) 

L6 = WF)icB(MT> 

L5 =L5+L6 
A(MT)= L6 
CNL = L5-CN 
GO TO 3 

780 WF=:( AG< IQ+1 )“AG( IQ-1 > )/2 

L6 = WFJicBdQ) 

L5 = L5+L6 
A(IQ)= L6 
3 CONTINUE 
890 C3 = 0 

DO 40 lY = IfMT 
IF (lY.GT.dHT+D) GO TO 980 
C5= (BdY)>KAH<IY)) 

AHdY)=C5 
C3 = C3+C5 

IF (IY*EQdHT) GO TO 970 
GO TO 40 . 

970 CCP = C3 
GO TO 40 

980 C5=(BdY>>}cAHdY)) 

AHdY) = C5 
C3 = C3 + C5 
40 CONTINUE 

PT-PPP 

IF(CLl.EQ^Od GO TO 196 
WRITE (7fl65) 

165 FORMAT(////dOX>' UPPER SURFACE") 

WRITE <7fl94) 

194 FORMAT (2X>" %CH0RD"f3Xf" CP LOCAL' r3Xf' CN LOCAL 'f 3X> 
C ' CC LOCAL' f3Xf' CM LOCAL') 

196 C9=0* 

C8=0* 

DO 5 S=lyMT 

IF (S,EQdIHT+l)) 60 TO 1240 
IF (S^GTdlHT+D) GO TO 1270 
BB(S) = (-A(G)>KAG(S)) + <AH(S))KAJ<S)) 

C9 = C9+BBCS) 

TIT = AG(S)5!aOO 
IF(CLl»EQ*Od GO TO 5 

WRITE(7d85) TIT^B(S) >A(S) >AH(S) »BB(S) 

GO TO 5 

1240 IF(CLl»EQ*Od GO TO 1270 
WRITE<7d95) 

195 F0RMAT<///10Xf ' LOWER SURFACE') 

WRITE (7f194) 

1270 BB(S) = (A(S))KAG<S) ) + (AH(S))KAJ(S) ) 

C8= C8-BB<S) 

TIT = AG(S)>SJ100 
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FORCE 


IF(CL1.EQ»0») GO TO 5 

WRITE (7^185) TITrB(S)»A(S)»AH(S)r-BB(S) 

185 F0RMAT<4X»F4»lr4F12»4) 

5 CONTINUE 

STl = C8+C9 

IF<CL1»E0,0) GO TO 209 
WRITE(7»205) 

205 F0RMAT(//18Xf 'PRESSURE' »2X» 'SUCTION' »2X» 'TOTAL V) 
WRITE(7»206)CNrCNLrL5 

206 F0RMAT<12X?'CN'f2X»3F9»4) 

WRITE (7f207)CCPfST?C3 

207 F0RMAT(12Xr 'CC' r2Xy3F9*4) 

WRITE(7f208)C9rC8»STl 

208 F0RMAT(12Xr'CM'r2Xr3F9*4) 

209 AT = AN :5c 0» 017453 
CSl = COS(AT) 

SN = SIN(AT) 

CL = (L5t>CSl)-(C3«SN) 

CD = (C3:{{CS1) + (L5>!cSN) 

WRITE (7y225) 

225 F0RMAT(22Xr' WING PERF0RMANCE'/18Xr 'CL ' » lOXf 'CD' > 10X» 'CM' ) 
WRITE <7»235) CL^CDrSTl 
235 FORMAT (10 Xt3F12. 4) 

CALL CL0SE(2) 

END 
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SUBROUTINE SET 
DIMENSION DA(IOO) 

C0MM0N/0NE/NJfMT?NRrCLlfBlFPRlrAKlfAK3fANfR3fPP2f ITRN 
C0MM0N/TU0/RN(30)fRS(30)fW(30)fX<30)fP(40)»Q(40)fRD<20) 
C0MM0N/THREE/D(60) f WTY(30) f UBY(30) r WL<30) rB<50) f PD<48 f4) 
COMMON/FOUR/E (30 ) fH(30)fY (30) fG (30) fWK30)FXI < 30) 
C0MM0N/FIUE/U(30) FV(30)FDS(60)fDET(60)FCS(20) 

COMMON/SEUEN/PC ( 50 ) 

COMMON/EIGHT/TAB f AMP f IFN f I AM 

DATA Y(1)fY(2)fG(1)fG(2)fY(23)fY(24)fG(23)fG(24)/0*f0*f0.f 
C O* fO# fO* fO# fO*/ 

CALL ASSIGN(3f'PAD.DAT'f0f'0LD'ff> 

DEFINE FILE 3 (50f512fUf INR) 

CALL ASSIGN(2f'NPL*DAT'fOf'OLD'ff) 

DEFINE FILE 2 (50f 128fUf INR) 

C 

C DATA OUTPUT 

C 

NJl = NJ-4 
IWT=1 

IF(CL1^EQ*0*) GO TO 12 
WRITE(7f65)' 

65 FORMAT(/" TRANSDUCER OUTPUT'/^ CP VALUES CHANNELS 2-4V 

C 7Xf' 1'f7Xf' 2'f7Xf' 3'f7Xf" 4') 

WRITE(7f70) (Nf(PD(NfJ)fJ=1f4)fN=1f48) 

70 F0RMAT(I3f4F944) 

URITE(7f41) ITRN 

41 FORMAT(//" RUN 'fI4f' OUTPUT V12Xf' EXT VEL* 'f7Xf' MOVEMENT' f 
C 8Xf' Y CO-ORD') 

DO 45 J == IfNJ 

IF(J^LT^3*0R*J*6Tt22) GO TO 63 
WMT = -WTY(J-2) 

WMB = -WBY(J-2) 

GO TO 62 
63 WMT = 0^0 

WMB = OtO 

62 WRITE (7f50) D( J) fE( J) fH( J) fY( J) fG( J) f WMT f WMB 
50 F0RMAT(7F8*4) 

45 CONTINUE 

WRITE(7f300)ITRN 

300 F0RMAT(//' RUN^fI4f' DATA'/ISXf' TOP WALL') 

WRITE(7f320) 

320 F0F^MAT(/5Xf ' JACK'fSXf' MACH NO*') 

TMN=0 

BMN=0 

CALL ASSIGN (6 F 'WALM* DAT' fOf 'OLD' f f ) 

DEFINE FILE 6 (20f256fUf INR) 

WRITE(7f21) 

21 FORMAT(//'$ MACH NO* RECORD (1-20) = ') 

READ(7f22)IMR 

22 F0RMAT(I6) 

DO 23 J=2f21 
X(J):=P(J) 

X(J+21)=Q(J) 
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C0NT.INUE 

WRITE(6'IMR) (X(J) r J=l»42) 

CALL CL0SE(6) 

GMN=0 

DO 170 J = 1»NJ1 

ISl = RS(J) + ( (Y( J+2) )*PC(J) ) 
DA(J)=IS1 
IP = RS(J) 

YM=P(J+1) 

WM = -WTY(J)-DET(J) 

IF (DA(16) ♦LT*123* ) DA(16)=123. 

IF (DA<17) ♦LT*178. ) DA<17)=178. 

IF (DA(18) .LT*161» ) DA(18)=161* 

IF (DA(19) ♦LT.140. ) DA(19)=140» 

IF (DA(20) .LT.216» > DA(20)=216* 
IF<CL1,EQ*0* )G0 TO 170 
URITE(7»840)JrUMyYMfY(J+2) rIPf ISl 
F0RMAT(7Xf I2»llXr3Fll«4,2(7XrI3)) 
IF< J»LT.5.0R. J»GT.13)G0 TO 170 
IBJ-J+20 

WRITE(7r810)P(IBJ) 

FORMAT <13X»F11 *4) 

CONTINUE 

IF(CLl.GT.O* ) URITE(7»20) 
F0RMAT(///15Xf ' BOTTOM WALL') 

DO 180 J = IrNJl 

ISl = RN(J) - ( (G( J+2) )>KPC( J+NJl) ) 

DA<J+NJ1)=IS1 

IP = RN<J) 

WM=-WBY ( J ) +DET ( J+N J1 ) 

NIJ=J+16 

YM=Q(J+1) 

IF <DA(36),LT.209*> DA(36)=209* 

IF (DA<37) »LT,143. ) DA(37)=143. 

IF (DA<38)J_T.136*) DA(38)=136» 

IF (DA<39) ♦LT*183* ) DA<39)=183» 

IF (DA<40) ♦LT.270* ) DA<40)=270» 
IF<CL1,EQ,0* )G0 TO 180 
WRI TE ( 7 » 840 ) J r WM ? YM » G ( J+2 ) » IP x I SI 
GO TO 180 

IF(J«GT*19) GO TO 180 
AMN=(P( J-fl)-Q(J+l) )/2 
GMN=GMN+(AMN^<AMN) 

TMN-=TMN+ < P < J+1 ) 5KP ( J+1 ) ) 

BMN=BMN+(Q( J+1)>KQ(J+1)) 
IF(J.LT*5.0R« J,GT»13) GO TO 180 
IBJ=J+20 

WRITE<7i-ai0)0<IBJ) 

CONTINUE 
KT = 25!<NJ1 
DO 190 J = 1»NJ 
DA(J+KT) = E(J) 

KS = KT 4 NJ 
DO 200 J = IrNJ 
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200 DA(J+KS) = H(J) 

IF(CL1*EQ*0* )WRITE(5>30) IFN 
WRITE<7f30) IFN 

30 FORMAK/' WALL & MODEL OUTPUT RECORD NO* = 'yI4//50<' 

35 F0RMAT(I2) 

ND = 2 * (NJ+NJl) 

WRITE (3 'IFN) <DA( J) »J=lfND> 

MT2 = MT + 2 
IF(IWT*EQ*1) MT2=MT 
WRITE(2'IFN) <B(J)»J=1»MT2) 

CALL CL0SE(2) 

CALL CL0SE(3) 

RETURN 

END 


')//) 
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SUBROUTINE WALL 
C 

C ONLINE WALL EXERCISER 
C 

COMMON/EIGHT/TABf AMPr IFNf lAM 
NWJ=40 

DIMENSION PC<50)»IM(50),ADC<50)rPV<50) 

DIMENSION PL(IO) 

DATA PL/130* r 178* f 161* f 140 ♦ f 216. f 209* f 143* r 136. f 183* f 270*/ 
CALL ASSIGN (3 > 'PAD. DAT' fOf 'OLD' f f ) 

DEFINE FILE 3 (50r512FUf INR) 

READ(3'IFN) (P0< J) f J=1 rNWJ) 

CALL CL0SE<3) 

C 

C CALCULATE MOVEMENT 

C 

IF=1 

IP=5 

CALL INIT 
50 IW=0 

ICNT=0 

IS=0 

IT0L=3 

DO 25 L=1»NWJ 
IAI=L+23 

IJP=CVSWG < I ADC ( lAI ) ) 

IPV=IJP5K* 24414 

IMOVE:=ADC<L)-IPV 

ADC(L)=IPV 

IF(IPV.GE.980) WRITE(5>200) L 
200 FORMAT</' POT LIMIT REACHED ON JACK 'fI4) 

IF<IPV.GE.980) STOP 
IFCL.LE.IS) GO TO 60 
IF(L.6T.20*AND.L.LE.35) GO TO 60 
IF(L.EQ.36) IW=15 
NDJ=L“(15+IW) 

IF<PV(L) .LT.PL(NDJ)) PV(L)=PL(NDJ) 

60 IPD=PV(L)-IP 

IMV=IPV-IPD 
IJ=L 

IF(IP*EQ*0.AND.IMV.GT.0)IMV=0 
IF(L.GT.20) IMV=:-IMV 
IF(L.GT.20) GO TO 90 
IEV=IJ/2 
IEV:=IEV+IEV 

IF(IEV.NE.IJ) IMV=-IMV 
90 IFdF.EQ.DGO TO 300 

IF(L.EQ.1.0R.L.EQ.21) GO TO 110 

IFaADS<IM(L>).GT.50.AND.IABS<IM0VE).LE.10> GO TO 400 
IF(IS.GE.39)ICNT==ICNT+1 . 

IF(IS.GE.39*AND.ICNT.EQ.6> GO TO 500 
110 IF<IABS<IMV)-IABS<IM(L)).GT.50) GO TO 500 
300 . IM(L)=IMV 

IF<IABS(IMV).LE.ITOL) GO TO 40 


Figure 4, 9.1 



0066 

0067 

0068 

0069 

0070 

0071 

0072 

0074 

0075 
0077 

0079 

0080 
0081 
0082 
0083 
0085 

0087 

0088 
0089 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
WALL 


0001 


0002 

0003 

0004 

0005 

0006 
0007 

0009 

0010 
0011 
0012 
INIT 


GO TO 35 
40 IS=IS+1 
IMV=0 

35 CALL MOUE(IJrlMO) 

25 CONTINUE 

IF=0 

IF(IS4EQ»40) go to 70 
CALL START(POflM) 

IF(IS*LT*40) GO TO 50 
70 IF(IPtEQ*0) GO TO 100 

IP=0 

GO TO 50 
100 IW=0 

DO 85 J=16>40 

IF( J*GT*20tAND*J*LTt36) GO TO 85 

IF(J*EQ*36)IW=15 

ILIM=J-(15+IW) 

IE=PV(J)-PL<ILIM) 

IF<IABS(IE).LE^ITOL) WRITE<5f80) J 
80 FORMAT(' LOWER POT LIMIT REACHED ON JACK ' flA) 

85 CONTINUE 

GO TO 250 

400 URITE(5rl50) L 

150 FORMATC//' JACK 'flAf' JAMMED!!!!!!!!!!!") 

GO TO 250 

500 WRITE<5>160) L 

160 FORMAT(//' WALL OUT OF CONTROL AT JACK " f I4/40< "*" ) ) 
250 RETURN 

END 


SUBROUTINE INIT 
C 

C INITIALISE TSWT CONTROL SYSTEM 

C 

NWJ=40 
IDI=128 
DO 10 J=lfNWJ 
ICOM=J+IDI-fO 

CALL DI0(IC0M>IDIf INPUT) 
IF<INPUT*NE»0) WRITE<5»30) J 
30 FORMAT(/' JACK 'rlAf' I/O ERROR') 
10 CONTINUE 

RETURN 
END 
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SUBROUTINE START(PVrlM) 

DIMENSION PV(50)fIM<50) 

ISUM=0 

DO 10 J=21f40 
ISUM=ISUM+IABS(IM(J)) 

IAM=ISUM/20 

RCT=(IAM/2)*25* 

IF(RCT*GT*6004) RCT=600- 
IF(RCT.LT*25*) RCT=25* 

IF<RCT.GT*50*) CALL DIO ( 129 f 128 f INPUT ) 
IF<RCT*GT.50») CALL DIO( 149f 128f INPUT) 

START ALL MOTORS 

IPSA=45 

IDI=0 

IC0M=IPSA+IDI+4096 
CALL DIO(ICOMflDIflNPUT) 

WAIT FOR WALL TO MOVE 

IEND=0 

CALL SETRCSrOfRCTflEND) 

CALL LWAIT<0rIEND) 

STOP ALL MOTORS 

CALL DI0<45fIDIf INPUT) 

RETURN 

END 
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SUBROUTINE MOVE(IJflMV) 

IF<IMV-0) ISylOfS 

IDR=3072 

GO TO 20 

IDR=0 

GO TO 20 

IDR==2560 

IDI=128 

ICOM=IJiIDI+IDR 

LOAD COMMAND 

CALL DIO<ICOM»IDIfINPUT) 

RETURN 

END 
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SUBROUTINE DIO ( ICOMr IDI , INPUT) 

CALL IPOKE< "167774 r 45) 

CALL IPOKE( -167774fIC0M) 

IF(IPEEK< -167770) ♦GE^O) GO TO 10 
CALL IPOKE( -167770^0) 

IF<IDI*EQ^O) GO TO 20 
IF<IPEEK(- 167770) *NE^ 128) GO TO 30 
INPUT=IPEEK< - 167772) 

CALL IPOKE< -167770f0) 

RETURN 

END 
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MAIN CONTROL PROGRAM 

COMMON/ONE/N J f MT r NR f CLl f B1 f PRl r AKl r AK3 r AN f R3 f PP2 f ITRN 
C0MM0N/TW0/RN(30)fRS(30)fW<30)fX<30)fP<40)fQ<40)fRD< 20) 
C0MM0N/THREE/ri(60) fWTY(30> fWBY(30> fWL(30) fB(50) fPD(48f4) 
COMMON/FOUR/E < 30 )fH< 30) fY (30) fG<30)fUI <30 )fXI (30) 
C0MM0N/FI0E/U(30)fV(30)fDS(60) fDET(60)fCS<20) 

COMMON/SEOEN/PC ( 50 ) 

COMMON/EIGHT/TABfAMPfIFNfIAM 
CLl = 500t 
NR = 8 
NJ = 24 
WRITE(7 f5) 

FORMAKIOXf' UNIVERSITY OF S0UTHAMPT0N'/5Xf 
C ' TRANSONIC SELF-STREAMLINING WIND TUNNEL V20Xf8< ) 
WRITE(7f10) 

F0RMAT(/17Xf ' DATA REANALYSISV17Xf 15(/3icO ) 

CALL IDATE(IfJfK) 

WRITE(7f15) JfIfK 
FORMAT ( 19X F 12 f2('-'f 12)) 

WRITE(7f40) 

FORMAT</'$ RUN NO* = ') 

READ(7f45) ITRN 
F0RMAT(I4) 

WRITE(7f110) 

F0RMAT(/'$ 

READ(7f45) IFN 
WRITE(7f50) 

F0RMAT(/'$ 

READ(7f60) an 
F0RMAT(F9*4) 

URITE(7f70) 

FORMAT(/'i 
READ(7 f45) MT 
URITE(7f80) 

FORMAT(/'$ INPUT 1 FOR AUTO IN-FILE SELECTION -') 

READ(7f45) IAM 

WRITE(7f85) 

FORMAT</' INPUT AMBIENT CONDITIONS V'lfc TEMP (DEG*C) = ') 

READ(7f25) tab 

F0RMAT(F8*4) 

WRITE(7f90) 

FORMAT('$ PRES* (CM HG) = ') 

READ(7f25) AMP 
CALL DATA 
CALL REDUCE 
CALL WAS 
CALL STAR 
CALL SUME 

IF(MT*EQ*0) GO TO lOO 
CALL FORCE 
CALL SET 
END 


FILE NO* = M 
MODEL ALPHA (DEG) = ') 

NO* OF MODEL TAPS ? ') 


Figure 4.10.1 



0001 


PROGRAM RUN 

0002 


DIMENSION TUN<5) 

0003 


CALL ASSIGN<2»'RUN.DAT'»0»"0LD'»») 

0004 


DEFINE FILE 2 (6r256»U» IJR) 

0005 


READ<2'l)(TUN(J)»J=lr5) 

0006 


WRITE<5»25) 

0007 

25 

FORMAT<'$ INPUT 1 FOR RECORD NO. INCREMENT 

0008 


READ(5y35)INC 

0009 

35 

FORMAT (16) 

0010 


IAD-CVSUG(IADC<21) ) 

0011 


TUN(3)=(IAD-2045)/8.15 

0012 


IF(INC.EQ.l) GO TO 20 

0014 


WRITE(5»5) 

0015 

5 

FORMAT (//' GOOD MORNING'/" INPUT PRES.') 

0016 


. READ(5»45)TUN(4) 

0017 

45 

FORMAT <F8. 4) 

0018 


WRITE(5fl5) 

0019 

15 

FORMAT(/' INPUT NEXT RUN NO. % FILE NO.') 

0020 


READ(5,10)IRfIF 

0021 

10 

FORMAT (21 6) 

0022 


TUN(1)=IR 

0023 


TUN(2)=IF' 

0024 


GO TO 30 

0025 

20 

TUN(2)=TUN(2)+1 

0026 

30 

WRITE(2'1) <TUN( J) r J=1^5) 

0027 


WRITE(5r55) (TUN(J) r J=l,5) 

0028 

55 

FORMAT(/' RUN DATA STORED '/5F8.2) 

0029 


CALL CL0SE(2) 

0030 


END 

RUN 




Figure 4.11.1 



Figure 5 A sample output of the control software OFLEX 



UNIVERSITY OF SOUTHAMPTON 
TRANSONIC SELF-STREAMLINING WIND TUNNEL 

/|V /fV /fV /|V 

AUTO MODE 
22 - 1 - 8 ?, 

RUN NO. = 407 

MODEL ALPHA <DEG) = 4.0 

AMBIENT CONDITIONS 

TEMP = 22.25 PRES<CM HG) = 76.65 

ITERATION RECORD NO. = 31 

MACH NO. = 0.5945 

REYNOLDS NO. - 0.118538E+07 

RECORD = 30 

WALL CP ERROR 

TOP - 0.0061 BOTTOM - '0.0055 

RESIDUALS = 0.0287 -0.0042 0.0014 

WALL a MODEL OUTPUT RECORD NO. = 31 


ITERATION RECORD NO. = 32 

MACH NO. = 0.5948 

REYNOLDS NO. - 0.118590E+07 

RECORD = 31 

WALL CP ERROR 

TOP - 0.0026 BOTTOM - 0.0038 

RESIDUALS = -0.0163 0.0051 0.0000 

WALL & MODEL OUTPUT RECORD NO . = 32 


Figure 5.1 



ITERATION RECORD NO. = 33 


HACH NO. = 0.5966 


REYNOLDS NO. - 0.118846E+07 

RECORD = 32 

WALL CP ERROR 

TOP - 0.0029 BOTTOM - 

RESIDUALS = -0.0066 0.0101 -0.0017 


# WALLS STREAMLINED ♦ 


WALL % MODEL OUTPUT RECORD NO. = 33 


0.0028 



Figure 6. A sample output of the 

re-analysis software QRLEX 



RUN ORLEX 


UNIVERSITY OF SOUTHAMPTON 
TRANSONIC SELF-STREAMLINING WIND TUNNEL 

^ \!/ nW vV nI/ nV vV 
/|\ 

DATA REANALYSIS' 

>K )i( ;K ^ >K>K ^ :K >joK >lc :|c 
11- 2-82 


RUN NO* = 389 

FILE NO* = 21 

MODEL ALPHA (DEG) = 3*0 

NO* OF MODEL TAPS ? 50 

INPUT 1 FOR AUTO IN-FILE SELECTION -1 

INPUT AMBIENT CONDITIONS 
TEMP <DEG*C) =22*0 
PRES* (CM HG) = 76*32 

DATA INPUT FILE = >KADC*DAT 
>ic * )i< * JioK )jC)io}c>K >!o!o!<>!< :jo!oK ** -t >K * ^ ^ >K 5i< ^ *!{< >K Jic 

ITERATION RECORD NO* = 21 

#)}:>!{#)!( >i( * jjoK >K >!o!< ^ >K ^ 5io!c>!<>!olo!< >!oioi«K ** >K »: >K **>}: >!< * 


MACH NO* = 0*8038 


REYNOLDS NO* - 0*143046E+07 

FILE NO* ? 

1 FOR M=<*725 

2 FOR *725<M>*825 

3 FOR M=>*825 ANS = 2 

WALL CONTOURS RECORD = 20 


WAS COMPUTING NOW ! ! 



DELTA STAR CALCS » 


TOP WALL 
TAP NO. 

DU/DX 

MACH NO. 

D* 

DD* 

1 

0.3870 

0.8238 

0.0155 

-0.0005 

2 

-3.0340 

0.8057 

0.0227 

-0.0004 

3 

0.6563 

0.7998 

0 . 0290 

-0.0000 

4 

0.1942 

0.8092 

0.0344 

-0.0009 

5 

1.9132 

0.8012 

0.0395 

-0.0012 

6 

14.1533 

0.8127 

0.0403 

-0.0042 

7 

32.5273 

0.8366 

0.0386 

-0.0076 

8 

54.5744 

0.8975 

0.0349 

-0.0130 

9 

55.2920 

0.9884 

0.0302 

-0.0189 

10 

-2.2115 

1.0633 

0.0294 

-0.0218 

11 

-53.4752 

0.9817 

0.0344 

-0.0189 

12 

-37.5962 

0.9032 

0.0424 

-0.0130 

13 

-17.6038 

0.3758 

0.0483 

-0.0088 

14 

-10.6314 

0.8436 

0.0571 

-0.0038 

15 

-4.2687 

0.8159 

0.0672 

0.0013 

16 

-0.3833 

0.8144 

0.0735 

0.0029 

17 

0.1915 

0.8130 

0.0781 

0.0033 

18 

0 . 5.746 

0.8159 

0.0823 

0.0029 

19 

0.1913 

0.8173 

0.0861 

0.0021 

20 

0.0000 

0.8173 

0.0907 

0.0021 


BOTTOM 

WALL 




TAP NO 

. DU/DX 

MACH NO. 

D>ic 

DD* 

1 

0.4216 

0.8077 

0.0155 

-0.0009 

9 

Am 

-1.1125 

0.8060 

0.0223 

-0.0004 

3 

-1.5025 

0.7981 

0.0290 

-0.0000 

4 

-3.5534 

0.7935 

0.0353 

0.0004 

5 

-7.0914 

0.7721 

0.0437 

0.0025 

6 

-12.0694 

0.7533 

0.0508 

0.0063 

7 

-3.1566 

0.7367 

0.0542 

0.0084 

8 

8.0421 

0.7459 

0.0550 

0.0075 

9 

7.1233 

0.7563 

0.0550 

0.0059 

10 

6.3412 

0.7634 

0.0554 

0.0046 

11 

-3.0033 

0.7714 

0.0567 

0.0042 

12 

-8.1537 

0.7563 

0.0596 

0.0058 

13 

0.6964 

0.7520 

0.0622 

0.0063 

14 

8.0705 

0.7742 

0.0630 

0.0025 

15 

3.0736 

0.8001 

0.0634 

-0.0017 

16 

1.9875 

0.8001 

0.0664 

-0.0037 

17 

1.6074 

0.8149 

0.0697 

-0.0055 

18 

-2.0260 

0.8121 

0.0748 

-0.0054 

19 

-1.6459 

0.7998 

0.0806 

-0.0051 

20 

0.0000 

0.7998 

0.0861 

-0.0046 


UNIT REYNOLDS NO* = 357276.9 D« FPG = 0.0088 


Piaure> . 7 



WALL CP ERROR 


TOP - 

0*0072 

BOTTOM - 

0*0064 

RESIDUAL ERRORS 



X 

U/UFS 

V/UFS 


-18*0000 

0*0009 

-0.0000 


-17*0000 

0*0009 

-0*0000 


-16*0000 

0*0008 

-0*0000 


-15*0000 

0*0007 

-0*0000 


-14*0000 

0*0006 

0*0000 


-13*0000 

0*0004 

-0*0000 


-12*0000 

0*0003 

-0*0000 


-11*0000 

0*0002 

-0.0001 


-10*0000 

0*0001 

-0*0001 


-9*0000 

-0*0000 

-0*0001 

i 

-8*0000 

-0*0002 

-0*0001 


-7*0000 

-0*0004 

-0.0001 


-6*0000 

-0*0005 

-0*0000 


-5*0000 

-0*0005 

-0*0000 


-4*0000 

-0*0002 

0*0000 


-3*0000 

0*0005 

-0*0000 


-2*0000 

0*0013 

-0*0000 


-1*0000 

0*0021 

-0*0001 


0*0000 

0*0024 

-0*0001 


1*0000 

0*0021 

-0.0001 


2*0000 

0*0014 

-0*0001 


3*0000 

0*0008 

-0*0001 


4^000 

^^^005 

-A 


5*0000 

0*0003 

-0*0001 


6*0000 

0*0002 

-0.0001 


7*0000 

0*0001 

-0*0001 


8*0000 

0*0000 

-0*0001 


9*0000 

-0*0001 

-0*0001 


10*0000 

-0*0001 

-0.0000 


11*0000 

-0*0002 

-0*0001 


12*0000 

-0*0001 

-0*0001 


13*0000 

-0*0000 

-0*0001 


14*0000 

0*0001 

-0.0000 


15*0000 

0*0001 

-0*0000 


16*0000 

0*0001 

-0*0001 


17*0000 

0*0001 

-0.0001 


18*0000 

0*0001 

-0.0001 


MODEL ERRORS 


CP 

-1*4979 

0*0017 

-0*0001 

-0*0034 

-0*7490 

0*0022 

-0.0001 

-0*0044 

0*0000 

0.0024 

-0*0001 

-0*0048 

0*7490 

0*0022 

-0.0002 

-0*0045 

1*4979 

0*0018 

-0*0001 

-0*0036 

2*2469 

0*0013 

-0.0001 

-0*0025 

2*9959 

0*0008 

-0.0001 

-0*0016 

3*7449 

0*0005 

-0*0001 

-0*0010 

4*4938 

0*0004 

-0*0001 

-0.0007 

EFFECT 


DELTA 1 

ALPHA ERROR => -0*0045 

DEGREES 

-0*0005 

INDUCED CAMBER =* 0*0007 DEGREES 

-0*0013 

VEL* ERROR 

CP =» -0.0048 



AUERAGE 

= -0.0029 


0.0048 


Figure 6.3. 


4 WALLS STREAMLINED * 
44444444444404444444444 



au sECTJoa AtfALysjs 

9510 • 

RUN NO* a 389 
ALPHA a 3*00 


UPPER SURFACE 



%CHORO 

CP LOCAL 

CN LOCAL 

CC LOCAL 

CM LOCAL 

, 

0*0 

1.0227 

-0.0038 

0.0079 

0.0000 


0.7 

-0.7242 

0.0036 

-0.0062 

-0.0001 


1.0 

-1.0257 

0.0046 

-0.0024 

-o.ooot 

* : 

1.6 

-1.2861 

0.0064 

-0.0027 

-0.0002 


2.0 

-1.2979 

0.0066 

-0.0019 

-0.0002 


2.6 

-1.3172 

0.0070 

-0.0018 

-0.0002 

{ 

3.1 

-1.3824 

0.0068 

-0.0015 

-0.0002 


3.6 

-1.3838 

0.0146 

-0.0026 

-0.0006 


5.2 

-1.3318 

0.0270 

-0.0041 

-0.0015 


7.7 

-1.2563 

0.0314 

-0.0038 

-0.0025 

.. • 

10.2 

-1.1720 

0.0439 

-0.0040 

-0.0046 

* i 

15.2 

-1.1809 

0.0591 

-0.0044 

-0.0091 

, ■ 1 
1 

20.2 

-1.1601 

0.0580 

-0.0035 

-0.0119 

• • ; 

25.2 

-1.1279 

0.0563 

-0.0028 

-0.0143 

1 

30.2 

-1.1249 

0.0844 

-0.0031 

-0.0256 

I 

40.2 

-1.0382 

0.0779 

-0.0021 

-0.0314 

< 

45.2 

-0.7661 

0.0383 

-0.0004 

-0.0173 

•• 1 

50.2 

-0.5044 

0.0189 

-0.0001 

-0.0095 

1 

1 

52.7 

-0.4348 

0.0109 

0.0000 

-0.0057 


55.2 

-0.3941 

0.0099 

0.0001 

-0.0054 

. 

57.7 

-0.3548 

0.0089 

0.0001 

-0.0051 


60.2 

-0.3606 

0.0090 

0.0002 

-0.0054 

. 

62.7 

-0.3504 

0.0087 

0.0003 

-0.0055 


65.2 

-0.2664 

0.0067 

0.0003 

-0.0043 


67.7 

-0.2778 

0.0069 

0.0003 

-0.0047 


70.2 

-0.2736 

0.0102 

0.0007 

-0.0071 


75.2 

-0.2361 

0.0118 

0.0010 

-0.0088 


80.2 

-0.1737 

0.0088 

0.0010 

-0.0070 


85.3 

-0.0976 

0.0049 

0.0007 

-0.0041 

. 

90.1 

-0.0143 

0.0007 

0.0001 

-0.0006 

; 

95.2 

0.0626 

-0.0031 

-0.0007 

0.0029 

* 


0.1394 

-0.0034 

0.0003 

0.0034 

1 

• 1 


LOWER SURFACE 




i 

3CH0RD 

CP LOCAL 

CN LOCAL 

CC LOCAL 

CM LOCAL 

-1 

0.0 

1.0684 

0.0053 

0.0098 

0.0000 

1.0 

0.8331 

0.0074 

0.0098 

-0.0003 

t 

1.8 

0.5949 

0.0127 

0.0059 

-0.0004 

• * .<j 

5.3 

0.1998 

0.0085 

0.0026 

-0.0005 


10.2 

-0.0488 

-0.0024 

-0.0005 

0.0003 

\ 

15.3 

-0.1493 

-0.0074 

-0.0009 

0.0012 


20.2 

-0.2291 

-Or. 0142 

-0.0009 

0.0029 


27.7 

-0.3163 

-0.0275 

-0.0000 

0.0076 

' ♦ 

37.6 

-0.3310 

-0.0288 

0.0017 

0.0107 


45.1 

-0.2158 

-0.0162 

0.0020 

0.0072 

^4 

52.6 

-0.0488 

-0.0037 

0.0006 

0.0019 


60.1 

0.0888 

0.0067 

-0.0011 

-0.0040 

-i 

■i 

• 1 

67.5 

0.3073 

0.0179 

-0.0030 

-0.0121 

71.7 

0.3629 

0.0318 

-0.0039 

-0.0228 

1 

V 

85.1 

0.4490 

0.0412 

-0.0036 

-0.0351 

• 1 

90.1 

0.4430 

0.0221 

-0.0000 

-0.0200 


95.1 

0.3825 

0.0190 

0.0005 

-0.0180 

j 


0.1311 

0.0032 

0.0003 

-0.0032 


PRESSURE SUCTION TOTAL 


I 


CN 

0.6319 

0.0755 

0.7075 

CC 

-0.0350 

0.0193 

-0.0157 

CM 

-0.1870 ■ 

-0.0844 

-0.2715 


WING 

PERFORMANCE 


CL 

CD 

CM 


0.7073 

0.0214 

-0.2715 


Figure 6.4. 



TRANSDUCER OUTPUT 


CP VALUES CHANNELS 2-4 



1 

o 

3 

4 

1 

26*5666 

-0*0463 

0*089.1 

-0*0107 

2 

3*6152 

-0*0063 

0*1375 

-0*0069 

3 

42*8183 

0*0038 

0*1788 

0*0075 

4 

49*5844 

-0*0171 

0*1154 

0*0176 

5 

5&*-2891 

(Vr^006 

^t33^1 

- 0 ,^649 

6 

55*5544 

-0*0248 

0*1037 

0*1066 

7 

26*5003 

-0*0777 

0*1504 

0*1431 

8 

55*9856 

-0*2113 

0*1702 

0*1230 

9 

57*4449 

-0*4048 

0*1594 

0*1015 

10 

57*4781 

-0*5595 

0*8655 

0*0858 

11 

56*3836 

-0*3903 

0*0904 

0*0680 

12 

54*6921 

-0*2219 

0*0778 

0*1014 

13 

26*5334 

-0*1621 

0*1490 

0*1108 

14 

52*8015 

-0*0915 

0*1700 

0*0617 

15 

53*0005 

-0*0229 

0*1365 

0*0120 

16 

52*5362 

-0*0197 

-0*2664 

0*0120 

17 

52*1382 

-0*0164 

-0*2778 

-0*0207 

18 

52*0719 

-0*0227 

-0*2736 

-0*0144 

19 

26*6993 

-0*0259 

-0*2361 

0*0125 

20 

50*1150 

0*0341 

-0*1737 

0*1953 

21 

43*9792 

1*0152 

-0*0976 

-0*0094 

22 

38*0755 

0*0032 

-0*0143 

0*0507 

23 

36*9146 

-0*0306 

1*1970 

0*1238 

24 

35*9860 

-0*1368 

0*1394 

0*1330 

25 

26*9978 

1*0057 

0*0876 

0*0514 

26 

35*0905 

-0*2630 

0*1088 

0*1102 

27 

35*2231 

-0*5160 

0*1518 

0*0619 

28 

34*9910 

-0*7459 

0*1370 

1*0440 

29 

2*5870 

-0*2625 

0*1574 

0*0932 

30 

7*8605 

-0*2378 

0*1383 

0*0945 

31 

26*5334 

-0*2352 

0*1550 

0*0680 

32 

13*2004 

-0*3198 

0*1407 

0*0623 

33 

22*0559 

-0*1557 

0*1401 

0*1152 

34 

27*6279 

-0*1036 

0*1227 

0*0365 

35 

29*9165 

-0*0356 

0*1471 

0*0107 

36 

31*7075 

-0*0254 

0*1513 

0*1195 

37 

26*5666 

-0*0146 

0*1070 

-0*0069 

38 

33*6643 

-0*0140 

0*3073 

-0*0075 

39 

33*9960 

-0*0317 

0*3629 

0*0000 

40 

31*4090 

0*0343 

0*4490 

-0*0289 

41 

27*6279 

-0*0153 

0*4430 

0*0964 

42 

24*5434 

-0*2015 

0*3825 

0*1291 

43 

26*5334 

-0*6707 

0*1311 

0*0768 

44 

11*3099 

-0*2734 

0*1413 

0*0611 

45 

1*6252 

0*0057 

0*1131 

0*0504 

46 

9*7842 

-0*0922 

0*0778 

0*1121 

47 

1*6252 

-0*0165 

0*0892 

0*0787 

48 

0*2653 

-0*0540 

-0*0742 

0*0964 



I 


RUN 389 OUTPUT 

EXT VEL » 

MOVEMENT 

Y CO-ORD 

0*0000 

0*0000 

0*0000 

0*0000 

0*0000 

0*0000 

0*0000 

3*0000 

0*0000 

0*0000 

0*0000 

0*0000 

0*0000 

0*0000 

6*0000 

0*0117 

0*0013 

0*0007 

- 0*0005 

0*0154 

0*0353 

9*0000 

0*0008 

0*0011 

0*0005 

- 0*0003 

0*0179 

0*0471 

12*0000 

- 0*0007 

- 0*0021 

- 0*0001 

0*0000 

0*0546 

0*0881 

15*0000 

0*0048 

- 0*0056 

- 0*0000 

0*0004 

0*0867 

0*1115 

18*0000 

0*0018 

- 0*0170 

- 0*0003 

0*0012 

0*1444 

6*1544 

20*0000 

0*0119 

- 0*0279 

- 0*0004 

0*0008 

0*2857 

0*1460 

21*0000 

0*0255 

- 0*0408 

0*0004 

- 0.0000 

0.3369 

0*1057 

22*0000 

0*0602 

- 0*0386 

0*0014 

- 0*0014 

0.4937 

0*0711 

23*0000 

0*1059 

- 0*0394 

0*0029 

- 0*0030 

0.5257 

- 0*0047 

24*0000 

0*1359 

- 0*0456 

0*0041 

- 0*0039 

0.5652 

- O . O 082 

25*0000 

0*1034 

- 0*0247 

0*0034 

- 0*0015 

0.5882 

- 0*0118 

26*0000 

0*0681 

- 0*0262 

0*0014 

0.0000 

0.4967 

- 0*0083 

27*0000 

0*0506 

- 0*0301 

0*0006 

0*0005 

0.4577 

0*0461 

29*0000 

0*0267 

- 0*0189 

0*0006 

0*0004 

0.3369 

0*0532 

32*0000 

0*0067 

- 0*0033 

0*0005 

0*0011 

0.2595 

0*0637 

35*0000 

0*0069 

- 0*0025 

0*0001 

0.0015 

0.2430 

0.0118 

38*0000 

0*0042 

0*0057 

0*0004 

0*0015 

0.2060 

0.0461 

41*0000 

0*0071 

0*0045 

0*0002 

0*0016 

0.1725 

0.0200 

44*0000 

0*0070 

- 0*0039 

0*0007 

0*0021 

0.1389 

- 0.0130 

47*0000 

0*0070 

- 0*0039 

0*0007 

0*0023 

0.0863 

0.0439 

50*0000 

0*0074 

- 0*0023 

0*0000 

0.0000 

0.0000 

0.0000 

53*0000 

0*0074 

- 0*0023 

0*0000 

0*0000 

0*0000 

0.0000 


Pimivci C. C.‘ 



RUN 389 DATA 

TOP WALL 
JACK MACH NO* 


MACH NO* RECORD (1-20) = 20 


1 


0*0158 

0*7634 

0*0007 

467 

467 

2 


0*0183 

0*7714 

0*0005 

566 

566 

3 


0*0546 

0*7563 

-0*0001 

551 

550 

4 


0*0875 

0*7520 

-0*0000 

611 

610 

5 

0*8101 

0*1456 

0*7742 

-0*0003 

639 

638 

6 

0*8256 

0*2899 

0*8001 

-0*0004 

692 

691 

7 

0*8748 

0*3445 

0*8001 

0*0004 

746 

746 

8 

0*9339 

0*5067 

0*8149 

0*0014 

890 

891 

9 

1*0569 

0*5445 

0*8121 

0*0029 

922 

924 

10 

1*1771 

0*5870 

0*7998 

0*0041 

941 

944 

11 

0*9315 

0*6070 

0*7998 

0*0034 

933 

935 

CM 

0*9200 

0*5097 

0*9032 

0*0014 

896 

897 

13 

0*9189 

0*4665 

0*8758 

0*0006 

869 

869 

14 


0*3407 

0*8436 

0*0006 

870 

870 

15 


0*2582 

0*8159 

0*0005 

775 

775 

16 


0*2401 

0*8144 

0*0001 

318 

318 

17 


0*2027 

0*8130 

0*0004 

362 

362 

18 


0*1696 

0*8159 

0*0002 

339 

339 

19 


0*1368 

0*8173 

0*0007 

300 

300 

20 


0*0842 

0*8173 

0*0007 

365 

365 


1 

BOTTOM WALL 

0*0345 

0*8077 

-0*0005 

583 

583 

2 

0*0467 

0*8060 

-0*0003 

502 

502 

3 

0*0881 

0*7981 

0*0000 

523 

522 

4 

0*1119 

0*7935 

0*0004 

531 

530 

5 

0*1569 

0*7721 

0*0012 

495 

493 

6 

0*1523 

0*7533 

0*0008 

382 

381 

7 

0*1141 

0*7367 

-0*0000 

474 

474 

8 

0*0786 

0*7459 

-0*0014 

509 

510 

9 

0*0012 

0*7563 

-0*0030 

568 

570 

10 

-0*0035 

0*7634 

-0*0039 

617 

620 

11 

-0*0076 

0*7714. 

-0*0015 

576 

577 

12 

-0*0025 

0*7563 

0*0000 

593 

592 

13 

0.0524 

0*7520 

0*0005 

543 

542 

14 

0*0557 

0*7742 

0*0004 

662 

661 

15 

0*0620 

0*8001 

0*0011 

649 

648 

16 

0*0080 

0*8001 

0*0015 

320 

318 

17 

0*0407 

0*8149 

0.0015 

242 

240 

18 

0*0146 

0*8121 

0*0016 

210 

208 

19 

0*0180 

^t7998 

-0> 002-1 

“2^1'1 

209 

20 

0*0393 

0*7998 

0*0023 

295 

293 


WALL & MODEL OUTPUT RECORD NO* = 21 


STOP “ 


Fiaure 6.7 



Wall data not taken at points 1,2,22-24 


Assumed origin of 



FIG. 7 SOFTWARE REPRESENTATION OF EACH FLEXIBLE WALL 



APPENDIX A. DIGITAL INPUT /OUTPUT PROTOCOL 


FOR TSWr CONTROL SYSTEM 

The transfer of digital information between the computer and wind tunnel 
hardware is an important part of the on-line control system. Digital I/O 
involves a complex interaction of system software and hardware. A command 
code has been devised to simplify the operation of the control system. The 
protocol of digital I/O are described in the following sections. 


Data Format 

■ Each packet of information sent to the wind tunnel consists of the data 
destination and operation plus the data itself. Information from the wind 
tunnel consists of data only. All sets of information are contained in a single 
16-bit word in the following manner:- 

Digital Output Word (Equivalent to a six digit octal number in binary code) 



Digital Input Word 


15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 



Device Addresses 


Each address consists of six bits of binary information which corresponds 
to a decimal number referred to as the software address. The following scheme 
has been chosen: 


Device 


Stepper Motors 


Pulse Sequence 
Generator 


Scanivalves (1) 
( 2 ) 

Encoders (1) 
( 2 ) 

Power Supplies 
etc. 


Software 

Address 

1 44 

45 


46 

48 

47 

49 

50 -»■ 63 


Binary Representation 
543210 543210 

OOOQOl-101100 

5 4 3 2 1 0 


543210 5 4 3210 


11 


111 


Data Operation 


The transfer of data takes three forms - read only, write only and write 
before read. These operations determine the type of data to be sent, if any. 
For example, the read only function requires no data from the computer other 
than the device address. 


Information on data operation is sent to the wind tunnel by adding a chosen 
software value to the decimal equivalent of the digital ouput word. This effecti- 
vely sets the required bits for correct information transfer. The software values 
are chosen thus: 


Data 

Information 


Binary Representation 


Read only 

Write only 

Write before 
read 


Digit 


15 8 7 6 

□ 033 

[0 I °l°l°l 



6 3 


Software Value 


-32768 + 128 = -32640 


0 + 0=0 
0 + 128 = 128 


- A2 - 




Data 


Actions at the wind tunnel are determined by the corresponding data 1 or 
2 bits of the output word. Data 3 bits of the input word allow checks to be 
made on system devices. 

The output data is transferred by the setting of data 1 and 2 bits, 
achieved by the modification of the digital output word as for the data oper- 
ating information. 

The software values were chosen thus: 


Function 

a) Motor direction 
Stop 

Forward-go 

Reverse-go 

Digit 


Binary 

Representation 
Data 1 
11 10 9 



Data 2 


b) Scanivalve Move 

-Home 

-Step on one 

Digit 

c) Pulse sequence generator 

-Start 

d) Motor Power Supply 

-Off 

-On 

Digit 


14 13 12 


0 

0 

n 


0 

0 

0 


3 

14 13 12 


0 

0 

14 13 12 

E 

0 

0 



0 

□ 


Software Value 


0 

2560 

3072 


20480 

16384 


4096 


0 

8192 
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Function 


Binary 

Representation 


Software Value 


e) 


Pulse sequence generator 
Increment step size 
(currently non opera- 
tional) . 


14 13 12 11 10 9 8 



No. 10s 100s 1000s 


0 -> 32512 


Input data will be in binary code to allow simple software manipulation 
The complete input word will always be read regardless of the quantity of 
information being transferred. 

The types of data are as follows: 

Device No. of bits of 

Information 


Scanivalve - encoder output 7 

- at home position 1 

Pulse sequence generator 

- finish pulse 1 

- step value 7 

System monitor 16 

Motor direction 3 


Command Coding 


Each command sent to the wind tunnel must be unambiguous and provide 
information on data destination and operation and the data itself. This is 
achieved by placing a decimal code number in binary on the 16 output lines. 
This command number N is then decoded by the wind tunnel hardware and some 
operation performed. 

The conmand number for each operation is determined by the code 


+ 


Data 

Value 



Device 


Data 

N = 

Address 

+ 

Operation 




Value 


- A4 - 






using software values only. For example, if the scanivalve (1) is required to 
step on one port, the following command number would be sent to the wind tunnel 


N = 46 + 0 + 16384 = 16430 


Summary 

The interface between wind tunnel hardware and the coii 5 >uter is by means of 
16 output lines and 16 input lines plus control lines. The voltage levels on 
the output lines are controlled by the software generated command numbers 
described above. The voltage levels on the input lines are controlled by software 
selected tunnel hardware. All transfer of data is accompanied by a handshaking 
procedure to ensure correct sequencing of communication operations. 
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